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SEISMIC  WAVES  GENERATED  BY  CHEMICAL  EXPLOSIONS 


Summary  and  Recommendations. 

This  report  describes  the  results  of  research  on  some 
fundamental  properties  of  seismic  waves  generated  by  chem¬ 
ical  explosions.  The  work  was  carried  out  between  August  1 
i960  and  July  31,  1963.  In  the  course  of  the  program,  seis 
mograms  from  some  160  small  charges,  one  500  pound  surface 
burst,  one  100  ton  surface  charge,  and  a  buried  charge  of 
10  tons  were  recorded  at  close  range.  The  experiments  in¬ 
volving  the  small  charges  were  designed  to  elucidate  a 
variety  of  aspects  of  the  generation  of  Beismic  waves  by 
explosive  sources.  The  principal  BubjectB  investigated 
were: 

1)  The  properties  of  source  functions  obtained  by 
phase  equalization  of  dispersed  surface  waves 
recorded  at  close  range; 

2)  The  effect  on  surface  waves  of  distributing  the 
charge  in  space  and  time; 

3)  The  effect  of  source  depth  on  the  properties  of 
the  corapreBsional  and  Rayleigh  waves; 

4)  The  effect  of  the  medium  at  the  shot  on  the  prop¬ 
erties  of  the  seismic  signal; 

5)  The  properties  of  the  horizontally  polarized  shear 
waves  generated  by  an  explosion. 
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Most  of  the  objectives  of  the  work  were  accomplished. 
In  some  cases  the  techniques  employed  were  not  totally  ade¬ 
quate  to  completely  solve  the  problem  under  investigation. 
On  the  other  hand,  some  interesting  areas  for  further  study 
are  suggested  by  some  of  the  results. 

Some  of  the  principal  conclusions  of  the  work  are: 

1)  Even  at  close  range,  the  spectrum  of  the  observed 
surface  waves  may  be  so  narrow  that  the  computed 
source  function  is  oscillatory.  The  peak  of  the 
spectrum  may  remain  at  the  same  frequency  over  a 
wide  range  of  charge  sizes,  suggesting  that  the 
shot-point  medium  and  the  geologic  layering  may 
exert  an  influence  on  the  radiated  surface  waves 
that  overrides  the  expected  effect  of  yield.  The 
difficulty  of  interpreting  the  calculated  source 
function  in  terms  of  simple  force  systems  acting 
at  the  origin  is  clearly  demonstrated  by  the  re¬ 
sults. 

2)  When  a  distributed  explosion  is  fired  so  that  there 
is  no  interaction  of  the  component  detonations  in 
the  non-linear  region  around  each  charge, simple 
superposition  will  give  the  radiation  of  the  Ray¬ 
leigh  wave.  Time  delays  associated  with  the 
spatial  distribution,  as  well  as  those  between 
detonations  must  be  taken  into  account.  The  ori¬ 
entation  of  maxima  and  minima  of  Rayleigh  wave 
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amplitudes  may  be  controlled  as  desired  by  this 
technique.  Superposition  is  simplified  in  this 
case  because  the  radiation  of  Rayleigh  ■waves  from 
a  single  point  charge  is  completely  symmetric  at 
the  test  site  at  which  the  experiments  were  con¬ 
ducted,  eliminating  the  necessity  of  considering 
an  azimuthally  dependent  amplitude  function  in  the 
motion  synthesis.  The  distribution  of  Love  waves 
around  a  point  charge  is  complicated,  so  that  nu¬ 
merical  superposition  of  the  output  from  a  number 
of  charges  is  difficult.  The  grosB  properties  of 
the  Love  waves  from  a  distributed  charge  seem  to 
be  explained  by  simple  superposition. 

3)  The  effect  of  shot  depth  on  the  body  wave  properties 
cannot  be  separated  easily  from  the  effect  of  chang¬ 
ing  material  properties  with  depth.  The  amplitude 
measured  directly  above  the  source  from  a  series  of 

0.5  pound  shots  in  clay  approximated  a  decay  law  of 
-1.5 

d  *  .  The  prolonged,  oscillating  P-motion  observed 
from  shots  at  various  depths  in  clay  is  the  result 
of  the  combined  effect  of  the  primary  reflections 
from  the  near-surface  boundaries  and  the  oscilla¬ 
tory  nature  of  the  P-wave  emanating  from  the  source. 

The  effeot  on  the  Rayleigh  wave  of  increasing  the 
shot  depth  is  to  decrease  the  amplitude.  High-gain 
filtered  records  from  shots  at  various  depths  in 
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clay  contain  long  period  motion  In  the  region  Im¬ 
mediately  over  the  source.  The  period  Is  that  of 
the  Rayleigh  wave  observed  at  distance.  The  pres¬ 
ence  of  long  period  motion  directly  over  the  source 
supports  a  hypothesis  similar  to  that  proposed  by 
Dlx  for  the  generation  of  long  waves  from  an  ex¬ 
plosion. 

4)  Horizontally  polarized  shear  waves  are  always  pro¬ 
duced  by  a  burled  explosion.  The  portion  of  the 
total  seismic  energy  that  Is  transmitted  In  this 
form  is  primarily  a  function  of  the  medium,  with 
shots  In  competent  brittle  materials  producing  pro¬ 
portionally  larger  shear  waves. 

Love  waves  are  radiated  with  marked  azimuthal 
variation  of  amplitude,  marked  by  nodal  lines  and 
clear-cut  reversal  In  polarity  of  the  entire  wave¬ 
form  In  different  directions  outward  from  the 
source.  The  waves  are  thought  to  be  produced  di¬ 
rectly  by  the  shot  as  a  result  of  non-linear  proc¬ 
esses,  such  as  cracking,  in  the  highly  stressed 
material  near  the  explosion. 

SH  motion  Is  also  produced  by  mode  conversion  at 
non-horlzontal  boundaries  and  apparent  SH  motion 
may  result  from  congressional  head  waves  refracted 
at  dipping  interfaces.  Short  period  transverse 
horizontal  motion  is  often  seen  beginning  with  very 


small  amplitude  at  the  same  time  as,  or  Immediately 
after,  the  onset  of  the  first  congressional  motion. 
The  origin  of  this  motion  is  unexplained.  It  is 
found  on  the  records  from  small  explosions  in  this 
research,  as  well  as  seismograms  of  events  from  the 
Nevada  Test  Site. 

The  shear  body  wave  propagating  directly  from  the 
source  has  been  recorded  at  locations  directly  above 
buried  shots.  It  is  identified  by  particle  motion 
and  velocity. 

Recommendat ions .  The  experiments  carried  out  under 
this  project  and  the  analysis  of  the  data  have  served  to  elu¬ 
cidate  a  number  of  properties  of  explosion-generated  seismic 
waves.  When  a  seismologist  has  obtained  a  record  of  the 
ground  motion  due  to  any  seismic  source,  he  knows  that  the 
record  contains  valuable  Information  concerning  the  physics 
of  the  wave-generating  mechanism,  and  wave  propagation,  as 
well  as  data  pertinent  to  the  determinatie  .  of  the  proper¬ 
ties  of  the  medium  through  which  the  waves  h  ve  traveled 
to  the  seismograph.  The  decision  as  to  the  best  manner  in 
which  to  proceed  to  extract  this  information  is  often  a 
difficult  one.  A  variety  of  methods  were  applied  in  this 
work,  and  the  results  are  of  value  in  evaluating  the  vari¬ 
ous  approaches. 

Even  though  the  records  have  been  analyzed  in  detail 
for  the  particular  phenomenon  under  study  in  each  case, they 
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undoubtedly  contain  much  more  Information  that  can  be  ex¬ 
ploited  In  future  work.  At  this  time,  there  seems  to  be 
no  point  In  further  extensive  experimental  programs  using 
small  charges  unless  a  true  breakthrough  in  instrumenta¬ 
tion  provides  a  means  of  gaining  new  data,  especially  in 
the  high-stress  region  near  the  shot.  On  the  other  hand, 
specific  experiments  designed  to  test  particular  hypotheses 
will  be  very  desirable,  and  these  small  charges  provide  an 
inexpensive  and  rapid  means  of  collecting  the  data.  Cur¬ 
rently,  under  Contract  AF  19(628)-l689,  model  studies  using 
very  small  explosive  charges  are  being  conducted,  and  these 
may  prove  to  be  a  superior  method  of  gathering  information. 

It  is  recommended  that  a  series  of  experiments  in 
which  the  charge  size  is  increased  several  orders  of  magni¬ 
tude  over  those  used  here  be  considered  under  the  VELA- 
Uniform  program.  In  particular,  a  sequence  of  one  thousand 
pound  spherical  or  short  cylindrical  high-explosive  charges, 
fired  at  depths  from  1000  feet  below  the  surface,  at  inter¬ 
vals,  up  to  a  shot  on  the  free  surface,  and  even  10  feet 
or  so  above  the  surface,  would  be  of  great  value.  The  site 
for  such  an  experimental  series  should  be  in  a  region  of 
flat  topography,  underlain  by  flat-lying  beds,  with  no  sig¬ 
nificant  structural  changes  over  an  area  of  many  square 
miles.  A  site  in  the  Great  Plains  seems  most  promising. 

If  such  a  series  is  conducted,  it  should  be  heavily  instru¬ 
mented,  and  will  require  the  cooperative  participation  of 
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several  laboratories. 

A  further  recommendation  is  that  all  underground  nu¬ 
clear  shots  be  recorded  with  three-component  seismographs 
at  short  range  and  many  azimuths.  The  value  of  such  data 
for  studying  wave  generation  is  made  abundantly  clear  by 
the  results  of  this  research. 
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CHAPTER  1 
INTRODUCTION 

Studies  of  fundamental  properties  of  seismic  waves 
generated  by  explosions,  and  of  some  of  the  parameters  af- 
feotlng  these  properties,  have  been  carried  out  by  the 
Department  of  Oeophysies  and  Oeophysloal  Engineering  of 
St.  Louis  University,  under  Contract  AP  19(604)-7402. 

An  understanding  of  the  mechanism  by  which  a  burled  ex¬ 
plosion  generates  seismic  waves,  »nd  the  manner  In  which  the 
variable  factors  Involved  In  the  detonation  affect  wave 
generation  1.  essential  to  the  solution  of  one  of  the  crit¬ 
ical  aspect,  of  the  VELA-Unlfora  effort,  vis.,  the  Identi¬ 
fication  of  a  recorded  seismic  event  a.  an  earthquake  or  an 
explosion.  In  full-scale  tests.  Instrumentation  to  measure 
the  response  of  the  earth  In  the  sone  of  elastic  behavior 
usually  must  be  so  far  out  that  ordinary  variations  In  geol¬ 
ogy  and  topography  between  the  source  and  the  recording  site 
Play  an  Important  role  In  determining  the  character  of  the 
seismogram.  In  the  faoe  of  a  number  of  undetermined  vari¬ 
able  parameters,  the  Isolation  of  the  effect  of  any  one  1, 
difficult.  While  small  chemical  explosion,  cm  never  be  a 
complete  substitute  for  full-scale  tests  In  seismologies 
research,  they  do  offer  the  possibility  of  scaling  down  the 
experiment  so  that  It  1.  practical  to  have  uniform  known 
geology  from  shot  point  to  seismograph. 

In  this  research,  over  a  period  of  three  years,  a 
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total  of  some  160  small  charges,  0.25  to  15  pounds  at  posi¬ 
tions  from  38  meters  below  the  surface  to  3  meters  above  it, 
have  been  fired  and  the  resulting  ground  motion  measured  at 
distances  of  0  to  250  meters.  All  seismograms  were  re¬ 
corded  on  the  earth’s  surface,  although  subsurface  measure¬ 
ments  would  have  been  desirable  in  some  cases.  The  media 
in  which  the  shots  were  fired  were  loess,  clay,  alluvium, 
sandstone,  and  limestone.  Details  of  the  experiments,  in¬ 
strumentation,  description  of  the  test  sites,  and  field  pro¬ 
cedures  are  presented  in  Chapter  II. 

One  of  the  criteria  suggested  for  distinguishing  be¬ 
tween  an  earthquake  and  an  explosion  is  the  character  of  the 
source  function  found  by  phase  equalization  of  dispersed 
surface  waves.  Although  details  of  the  time  history  of  the 
force  applied  at  the  source  cannot  be  recovered  in  this  way, 
the  initial  polarity  and  type  of  force  (impulse,  step,  etc.,) 
can  in  principle  be  determined.  The  properties  of  the 
source  functions  obtained  under  a  variety  of  conditions  for 
shots  in  soil  have  been  investigated,  and  are  described  in 
Chapter  III. 

In  an  earthquake,  especially  one  of  large  magnitude 
in  which  a  fault  length  of  hundreds  of  kilometers  may  be 
involved,  a  point  source  is  an  inadequate  representation 
of  the  vjave -generating  mechanism,  and  the  effect  on  the 
seismic  signal  of  the  moving  source  must  be  taken  into  ac¬ 
count  in  a  careful  analysis.  To  simulate  a  moving  source. 
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several  shots  were  fired  in  which  the  charge  was  distribu¬ 
ted  In  a  number  of  holes  and  fired  sequentially.  The  re¬ 
sulting  radiation  patterns  of  Rayleigh  waves  are  described 
In  Chapter  IV,  while  the  effect  on  Love  waves  is  presented 
In  Chapter  V  on  explosion-generated  SH  motion. 

Early  In  the  research.  It  was  established  that  horizon¬ 
tally  polarized  shear  waves  (SH  waves)  are  a  normal  result 
of  a  burled  explosion  fired  under  a  wide  variety  of  condi¬ 
tions.  While  this  result  Is  unfortunate  from  the  viewpoint 
of  distinguishing  between  explosions  and  earthquakes.  It  has 
stimulated  further  thought  by  several  Investigators  on  the 
nature  of  a  contained  explosion  as  a  seismic  source.  An  ex¬ 
plosion  Is  essentially  a  center  of  compression  with  spher¬ 
ical  symmetry,  but  this  model  Is  Inadequate  to  explain  the 
details  of  the  resulting  seismic  waves. 

The  effect  of  source  depth  on  the  generation  of  comp- 
presslonal  and  Rayleigh  waves,  and  some  observations  on  the 
effect  of  the  medium  around  the  shot  are  discussed  In  Chap¬ 
ter  VI. 

Two  large  chemical  explosions,  one  of  100  tons  of  TOT 
fired  as  an  experiment  of  the  Canadian  Defense  Research 
Board  at  the  Suffleld  Experimental  Station,  and  one  of  10 
tons  of  Tovex  fired  near  Poplar  Bluff,  Missouri,  were  re¬ 
corded  as  part  of  the  research  program.  The  results  are 
discussed  In  the  Appendix. 

The  assistance  of  the  staff  of  the  Suffleld  Experimental 


Station,  Ralston,  Alberta,  and  in  particular  Mr.  G.H.S. 
Jones,  seismologist,  in  carrying  out  the  observations  at 
their  test  site  is  gratefully  aclcnowledged.  Gratitude  is 
also  expressed  to  the  managements  of  the  Alpha  Portland 
Cement  Company  quarry,  Lemay,  lissouri,  and  Pioneer  Silica 
Company,  Augusta,  Missouri,  for  permission  to  carry  out 
experiments  in  their  quarries. 


CHAPTER  II 


THE  EXPERIMENTAL  PROGRAM 

All  of  the  field  experiments  carried  out  under  the  re¬ 
search  program  are  summarised  In  Table  1.  The  experiments 
are  organised  according  to  the  test  site  at  which  they  were 
conducted.  For  all  exoept  the  tests  at  Suf field  Experi¬ 
mental  station  (SES)  the  explosive  used  w.s  forty  per  cent 
gelatin  dynamite.  In  one-half  pound  sticks,  eight  Inches 
long  by  one  and  one-eighth  Inch  diameter.  At  SES,  Nitrone 
was  used,  in  the  Table,  the  column  headed  "Array"  describes 
the  disposition  of  the  Instrumentation  as  follows:  "P"  ln_ 
dlcates  a  profile  In  u„e  with  the  shot  point,  and  the  num¬ 
bers  following  give  the  range  of  shot-to-detector  distances 
severed  by  the  profile;  "A"  Indicates  azimuthal  coverage  at 
a  fixed  distance  from  the  shot  point,  and  the  numbers  fol¬ 
lowing  give  the  angle  subtended  by  the  array  and  the  distance 
Description  of  the  Test  SlteH. 

Four  test  sites  were  used.  A  wider  variety  of  media 
was  originally  planned,  but  detailed  knowledge  of  the  geol¬ 
ogy  and  the  properties  of  the  waves  at  a  site  are  of  such 
Importance  In  Interpreting  the  experiments,  it  was  Judged 
advisable  to  concentrate  the  effort  Into  these  few  sites. 

The  materials  selected  are  representative  of  media  In  which 
test  explosions  are  likely  to  be  fired,  although  data  for 
Shots  In  crystalline  Igneous  rocks  would  have  been  useful. 

The  results  for  the  shots  m  compact  crystalline  limestone 
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are  probably  applicable  to  chemical  explosions  In  Igneous 
intruslves  of  fine  to  medium  texture.  Because  the  equation 
of  state  of  the  material  In  the  range  of  hydrodynamic  be¬ 
havior  Is  important  in  determining  the  response  to  a  nu¬ 
clear  shot,  caution  must  be  observed  in  carrying  over  the 
results  for  carbonate  rocks  to  aggregates  of  feldspar  and 
quartz. 

Florissant  Test  Site.  The  Florissant  Test  Site  is  an 
open  field  on  the  grounds  of  St.  Stanislaus  Seminary,  Floris 
sant,  Missouri.  Located  in  T  47  N,  R  6  E,  in  St.  Louis 
County,  the  site  is  roughly  a  mile  from  the  Missouri  River. 
The  surface  material  is  loess,  underlain  by  Pleistocene 
silts  and  clays.  The  bedrock  is  Mississippian  limestone, 
at  a  depth  of  33  to  45  meters.  Clay  densities  of  1.2  to 
2.0  gm/cm^  occur  above  the  water  table,  and  about  2.2  gm/cm^ 
below.  At  about  25  meters  depth  the  clay  changes  from 
stiff  to  very  stiff  (Kisslinger,  1959). 

Alpha  Test  Site.  The  Alpha  Test  Site  is  the  quarry 
operated  by  the  Alpha  Portland  Cement  Company  in  Lemay, 
Missouri.  It  is  located  in  the  northeast  quarter  of  Sec¬ 
tion  26,  T  44  N,  R  6  E,  St.  Louis  County.  The  exposed 
limestone  is  the  St.  Louis  formation,  of  Middle  Mississip¬ 
pian  age.  The  rock  is  a  dense,  hard  limestone,  varying 
from  lithographic  to  finely  crystalline  in  texture.  Chert 
nodules  are  widespread.  Over  part  of  the  quarry  floor  the 
surface  rock  is  almost  purely  calcitic,  in  other  parts  it 
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is  dolomitio.  The  average  density  is  2.67  gnv/cm3.  (Klss- 
linger,  1953)* 

Av.guata  Test  Site,  The  Augusta  Test  Site  is  a  quarry 
in  the  St.  Peter  sandstone,  located  along  the  bluffs  of  the 
Missouri  River  in  T  44  N,  R  1  E,  St.  Charles  County,  Mis¬ 
souri.  The  St.  Peter  sandstone  is  a  very  pure  silica  sand 
of  Middle  Ordovician  age.  The  well-rounded  grains  are  poor¬ 
ly  cemented .  The  grain  size  is  uniform,  0.12  -  0.25  mm,  The 
density  is  1.8  to  2.0  gn^/cmS. 

In  order  to  avoid  interference  with  quarrying  opera¬ 
tions,  all  of  the  experiments  were  confined  to  one  bench, 
some  300  meters  long  and  20  to  50  meters  wide.  About  25 
meters  of  sandstone  lay  under  the  level  at  which  the  tests 
were  conducted. 

Suffield  Test  Site.  The  Suffield  Test  Site  is  on  and 
adjacent  to  the  Watching  Hill  Blast  Range  of  the  Suffield 
Experimental  Station  of  the  Canadian  Defense  Research  Board, 
Ralston,  Alberta.  It  is  located  approximately  25  miles 
northwest  of  Medicine  Hat,  in  T  17  NandT  4.8  N,  R  5  W,  At 
the  site,  approximately  64  meters  of  silts  and  clays,  inter- 
bedded  with  gravels,  overlay  the  Cretaceous  Foremost  forma¬ 
tion.  The  Foremost  is  made  up  of  brackish-water  deposits 
of  interbedded  light  gray  sandstones,  dark  gray  shales  and 
coal  seams,  pale  yellow  silts,  and  rusty-colored  concretion¬ 
ary  banks.  The  densities  of  the  Burface  silts  and  clays  are 
low,  less  than  2  gm/em3-  <Ru8#e11  Landes,  1940;  Crock- 
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ford,  19^9). 

Instrumentation. 

The  project  can  be  divided  into  two  ptoses  on  the  basis 
of  the  instruments  employed.  During  the  first  year,  five 
Sprengnether  three-component  portable  seismographs  were 
available.  For  the  last  two  years  (starting  with  shot  point 
48),  a  sixth  standard  Sprengriether,  a  very  low-magnification 
Sprengnether,  and  a  12  channel  low  frequency  electronic  re¬ 
fraction  system  were  added. 

The  Sprenr.nether  Three-Component  Portable  Seismograph. 
This  instrument  was  the  primary  tool  of  this  research  pro¬ 
gram.  It  is  a  simple  mechanical-optical  seismograph,  in¬ 
corporating  all  three  components,  the  timer,  and  camera  into 
one  easily  emplaced  box.  (Figure  l).  Some  modification  of 
the  standard  instrument  was  required  to  meet  the  needs  of 
the  project.  A  system  was  developed  so  that  all  units  could 
be  turned  on  from  a  single  remote  control  point  and  various 
methods  of  putting  the  shot  firing  time  on  each  record  were 
tried.  The  zero-time  device  waB  never  completely  reliable, 
and  it  is  fortunate  that  much  of  this  work  required  wave¬ 
forms  rather  than  very  precise  travel  times.  During  the 
last  two  years  of  the  work,  better  timing  was  obtained  by 
tying  the  first  breaks  to  the  records  from  the  high-gain 
refraction  system. 

The  natural  undamped  period  of  all  seismometer  elements 
is  0,75  second,  and  all  components  are  damped  to  0,55  crit- 
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FIGURE  I 

SPRENGNETHER  THREE-COMPONENT 
PORTABLE  SEISMOGRAPH 
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lcal  by  a  conventional  damping  vane  -  magnet  arrangement. 

All  but  the  very  low  magnification  instrument  are  equipped 
with  dual  magnification,  in  which ,by  simple  optical  doub¬ 
ling,  a  high  gain,  twice  the  low  value,  is  available. 

The  low  magnification  instrument,  in  which  the  mirror 
is  mounted  directly  on  the  pendulum  for  each  component,  has 
a  static  magnification  of  7.5.  This  proved  very  useful  for 
measurements  directly  over  shots  burled  as  shallow  as  three 
meters,  and  for  observations  close  to  the  100  ton  and  10 
ton  shots. 

The  standard  instruments  are  actually  of  two  different 
generations.  Three  of  them.  Numbers  3,  4  and  5,  are  old  in¬ 
struments  that  were  reworked  for  the  purpose  of  this  pro¬ 
gram.  They  are  characterized  by  a  higher  magnification  of 
the  vertical  component  than  the  horizontals.  The  vertical 
components  on  two  of  these  were  marked  by  high  frequency 
noise  which  was  eventually  traced  to  the  camera  motor.  This 
noise  is  on  the  records  for  almost  all  or  the  shots. 

When  the  case  within  which  the  seismograph  is  mounted 
is  subjected  to  strong  shock,  as  in  some  of  these  close-in 
experiments,  the  panels  of  the  case  vibrate,  and  these  para¬ 
sitic  vibrations  show  up  on  the  seismograms.  This  situation 
was  improved  in  two  of  the  instruments  by  lining  the  entire 
inside  of  the  case  with  Y-9052  Scotch  Foam  Vibration  Danc¬ 
ing  Tape.  A  number  of  reproducibility  tests  were  carried 
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out  (see  Table  I)  which  aided  in  distinguishing  between 
true  ground  motion  and  instrument  noise. 

The  static  magnifications  (Low)  of  all  components  in 
the  new  instruments  (Numbers  1,  2,  and  6)  are  nominally  50. 

More  precise  values  were  determined  by  tilt-test  calibra¬ 
tion.  For  the  older  instruments,  the  horizontal  components 
have  magnifications  of  about  90,  and  the  verticals  about 
145.  High  magnifications  are  twice  these  values.  When  the 
instruments  are  operated  on  high  magnification,  the  direc¬ 
tion  of  motion  on  the  vertical  congjonent  only  is  reversed 
from  the  sense  on  low  magnification. 

The  SIE  GTR  200  Low  Frequency  Refraction  System.  A  12 
channel  unit  Identical  to  that  developed  for  the  crustal 
studies  program  of  the  United  States  Geological  Survey  was 
purchased  from  Dresser  Electronics,  SIE  Division  (Warrick, 
et  al,  1961).  The  geophoneB  selected  were  1.0  cpB  units 
(EV-12)  built  by  Electro-Technical  Laboratories.  The  per¬ 
formance  of  the  system  was  satisfactory  for  the  close-in 
measurements  program. 

Six  horizontal  component  1  cps  geophones  (EV-17H)  from 
Electro-Technical  Laboratories  were  acquired  toward  the  end 
of  the  contract.  Because  of  a  failure  in  the  GTR  200  sys¬ 
tem  it  waB  not  possible  to  make  any  records  with  these 
seismometers.  They  will  prove  useful  in  future  observations 
to  be  made  under  a  continuation  of  this  research  under  an¬ 
other  contract. 
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The  response  characteristics  of  the  GTR  200  are  given 
In  the  paper  by  Warrick,  at  al,  cited  above.  Galvanometrlc 
recording  only  was  used. 

Data  Processing. 

The  field  records  were  converted  to  digital  form  almost 
as  a  routine  practice.  The  digitalization  was  accomplished 
with  a  system  developed  by  Data  Instruments  Co.,  of  North 
Hollywood,  California.  The  heart  of  the  system  1b  a  unit 
known  as  the  Telecordex  Type  8G,  which  accepts  the  output 
from  magnetic  shaft  position  encoders,  and  puts  out  digital 
information  to  a  typewriter  or  an  IBM  Key  Punch.  Two  scan¬ 
ning  devices  were  used,  one  a  micrometer  microscope  on  which 
shaft  position  encoders  were  mounted,  the  other  a  Type  099D4 
Data  Reducer  from  Data  Instruments.  The  Data  Reducer  has  a 
reading  precision  of  157  counts  per  centimeter  of  record. 

The  microscope  offers  a  precision  of  15700  counts  per  centi¬ 
meter  of  record,  and  was  used  for  most  of  the  work  on  the 
records  from  small  explosions.  The  microscope  also  offered 
optical  magnification  of  ten  times,  which  was  of  great  value 
in  interpreting  the  small  signals  in  many  of  the  tests. 

The  digitalization  process  converts  times  and  ampli¬ 
tudes  into  counts.  A  computer  program  was  developed  for 
the  IBM  1620  to  convert  the  counts  back  into  time  and  to 
ground  amplitudes.  The  amplitudes  were  computed  by  assum¬ 
ing  that  the  seismographs  were  responding  as  true  displace¬ 
ment  meters  for  the  frequencies  involved,  and  the  trace 
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amplitudes  were  divided  by  the  static  magnification.  Time 
computations  presented  a  special  problem  since  the  chart 
speeds  were  slightly  different  for  the  various  seismographs, 
and  were  somewhat  variable  for  any  one  instrument.  The 
positions  of  the  timing  lines  over  the  portion  of  the  record 
to  be  processed  were  read  into  the  computer,  and  the  program 
interpolated  the  position  of  the  point  on  the  trace  it  which 
a  reading  was  taken  between  two  adjacent  timing  lines,  and 
expressed  this  position  as  the  time  from  the  selected  zero. 

Those  records  that  were  to  be  subjected  to  Fourier 
analysis  were  digitalized  at  equal  time  intervals,  spaced 
closely  enough  that  the  folding  frequency  was  twice  the 
highest  frequency  of  interest.  All  that  was  required  from 
other  records  was  the  values  of  peak  amplitudes,  time  of 
arrival  of  selected  phases,  or  sufficient  detail  to  allow 
the  construction  of  particle  motion  diagrams. 

The  following  programs  for  the  IBM  1620  were  developed 
as  part  of  the  research: 

Ground  motion  from  digitalized  seismograms 
Plot  ground  motion 
Fourier  transform 

Fourier  synthesis  of  phase  equalized  ground  motion 
Phase  velocities  from  Fourier  transforms 
Impulse  response  for  given  dispersion  curve 
Source  functions  by  convolution 
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Synthesis  of  ripple-fired  shots 
Travel  times  for  all  body  wave  arrivals  in  three¬ 
layered  medium. 

Detailed  discussions  of  many  of  the  seismograms  re¬ 
sulting  from  this  series  of  test3  have  been  presented  in 
the  series  of  Semi-Annual  Technical  Summaries  that  have 
been  issued  during  the  course  of  the  project,  as  well  as 
in  scientific  papers  and  special  reports.  It  is  not  prac¬ 
tical  to  present  a  complete  review  of  all  experimental  re¬ 
sults  in  this  report.  Rather,  the  main  results  will  be 
presented,  along  with  sufficient  data  to  support  the  con¬ 
clusions,  and  references  will  be  made  to  previously  publish¬ 
ed  material  where  appropriate. 


SUMMARY  OF  FIELD  PROGRAM 
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CHAPTER  III 


SOURCE  FUNCTIONS  AND  SYNTHESIS  OF  MOTION 
BY  PHASE  EQUALIZATION 

The  mechanism  of  seismic  wave  generation  is  one  of  the 
fundamental  problems  of  seismology.  One  objective  of  this 
program  of  research  was  to  gain  a  clearer  understanding  of 
the  process  by  which  surface  waves  are  generated  by  explo¬ 
sions.  Except  in  the  most  simple  case  of  an  elastic  half¬ 
space,  the  notion  in  the  source  region  which  gives  rise  to 
the  surface  waves  observed  at  distance  is  distorted  by  dis¬ 
persion  during  propagation.  In  principle,  the  effects  of 
dispersion  and  other  propagation  effects  can  be  removed,  and, 
subject  to  a  number  of  assumptions,  the  source  motion  can 
be  recovered. 

Two  equivalent  phase  equalization  methods  have  been 
used  with  equal  success.  These  are  described  by  Sato  (1955* 
1956)  and  Aki  (1960b).  The  analyses  of  source  motion  from 
the  fundamental  Rayleigh  mode,  the  first  Rayleigh  shear  mode, 
and  Love  waves  have  been  carried  out. 

The  first  method  is  based  upon  the  simple  principle 
that  the  amplitude  and  phase  angle  of  the  Fourier  trans¬ 
form  of  the  trace  of  the  motion  at  any  observing  station 
is  the  result  of  the  modification  of  the  initial  amplitude 
and  phase  spectra  near  the  origin  as  a  result  of  phase 
shifts  and  amplitude  changes  caused  by  dispersion,  attenua¬ 
tion,  and  recording.  The  Fourier  amplitude  and  phase  spec¬ 
tra  of  the  observed  motion  are  the  basic  data  in  this 
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technique .  From  the  spectral  data  along  an  observed  pro¬ 
file  the  phase  velocity  and  Initial  phases  of  all  components 
are  calculated,  absorption  coefficients  determined,  and 
t  inally  the  notion  at  any  distance  calculated  by  resynthe¬ 
sis  of  the  Fourier  conponents. 

Assuming  cylindrical  spreading  of  the  wavefront,  the 
motion  at  any  distance,  r.  Is  given  by 

(1)  f(r,t)  =  I  2  j  A(r,u>)  F 

V  77  r  J0  008  (  <P  (^)-cot)dcu 

where:  r«  ■  distance  to  recording  station 

)  =  observed  Fourier  spectral  amplitude 
V  ( )  =  initial  phase,  and 

A(r,  ev)  «  attenuation  function. 

The  attenuation  function  represents  the  absorption 
losses  within  the  medium.  It  is  assumed  that  these  can  be 
expressed  by: 

(2)  A  (r,  oj  )  =  exp  (-R<*  ) 

in  which  c*  can  be  written  as  Btu",  where  B  is  an  absorption 
coefficient  and  n  defines  the  dependence  of  attenuation  on 
frequency. 

One  of  the  principal  assumptions  in  the  techniques  used 
for  calculation  of  the  source  functions  is  that  the  spectrum 
at  any  observation  point  expresses  the  spectrum  at  the 
source.  If,  however,  higher  frequencies  have  been  attenuat¬ 
ed  to  the  point  that  they  are  not  recorded,  no  absorption 
law  will  replace  then.  The  loss  of  the  high  frequency  en¬ 
ergy  will  make  it  impossible  to  arrive  at  a  short  rise  time 
for  the  computed  source  function,  and  if  the  observed 
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spectrum  Is  severely  band-limited,  the  calculated  motion 
must  be  oscillatory. 

In  the  second  method,  first  exploited  by  Aki,  the  im¬ 
pulse  response  of  the  layered  medium  is  cross-correlated 
with  the  field  record  to  obtain  the  motion  at  the  source 
(source  function).  If  the  phase  velocity  curve  is  known, 
a  theoretical  seismogram  at  any  distance  can  be  calculated 


(3)  g(t)  =  j  G(w)  exp  (iojt)  6oj 

where  G(6<j)  represents  the  combined  frequency  response  of 
the  earth  and  recording  system.  |o(a;)|lB  assumed  to  be 
unity  only  over  a  finite  frequency  band,  sayu^  to  oj  2, 
and  the  phase  of  0(6j)  includes  shifts  due  to  dispersion 
and  recording.  The  cross-correlation  of  a  theoretical 
seismogram,  g(t),  along  a  particular  wave  path  with  the 
field  record,  x(tj,  for  the  same  path  will  yield  the  source 
function  y  (nr ) ,  i ,e . 

OO 

(4)  y  (t)  =  fg(t)  x  (t+T  )dt. 

Jo 

It  has  been  shown  (Jones  and  Morrison,  1954)  that  this  cross- 
correlation  equalizes  the  phase  angle  of  the  record  x  (t) 
to  the  source  in  the  frequency  range  from  ^  ^  to  ^2* 

Test  Sites  and  Data. 


The  surface  wave  data  used  for  the  calculations  of  the 
source  functions  were  obtained  at  two  test  sites*  the  Flor¬ 
issant  site,  which  consists  of  approximately  120  feet  of 
clay-loesB  overlying  Mississippian  limestone;  and  the  test 
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site  at  the  Suffield  Experimental  Station,  which  has  about 
210  feet  of  clay  with  sand  and  gravel  lenses  overlying  the 
Cretaceous  foremost  sandy  shale  formation  (see  Chapter  II). 

Examples  of  the  recorded  ground  motions  which  were 
used  for  source  function  calculations  are  shown  in  Figures 
2  and  3»  Figure  2  is  a  profile  of  the  three  components  of 
motion  of  the  surface  waves  at  Florissant  (S.P.  6,  7).  The 
motion  was  recorded  on  Sprengnether  portable  three- 
component  blast  seismographs.  The  instruments  were  care¬ 
fully  aligned  with  the  shot  point  so  that  one  horizontal 
component  responded  to  radial  motion  and  the  other  to  tan¬ 
gential  motion.  The  vertical  and  radial  components  are 
dominated  by  the  fundamental  Rayleigh  mode  which  is  pre¬ 
ceded  by  prograde  elliptic  motion  (Kisslinger,  1959;  Gilbert, 
et  al  1962).  This  prograde  motion  is  thought  to  be  a  com¬ 
bination  of  7  (a  diffracted  congressional  wave)  and  the 
first  Rayleigh  shear  mode.  Both  the  prograde  and  retro¬ 
grade  events  show  normal  dispersion  across  the  profile. 

The  transverse  motion  is  identified  as  a  Love  wave,  and 
also  is  normally  dispersed. 

Figure  3  shows  the  recorded  motion  along  a  profile  at 
the  Suffield  test  site,  produced  by  a  10  lb.  charge  at  a 
depth  of  5.5  meters  (S.P.  44a).  The  vertical  and  radial 
component  motion  is  dominated  by  the  fundamental  Rayleigh 
mode  at  the  closer  distances.  This  event  is  preceded  by  a 
prograde  elliptic  event  which  has  been  identified  as  the 
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firat  Rayleigh  shear  mode  M21.  This  identification  ia 
based  upon  the  particle  motions.  Figure  4,  and  the  phase 
velocity  curve.  Figure  8.  The  development  of  this  event  as 
it  separates  from  the  fundamental  Rayleigh  wave  with  increas¬ 
ing  distance  is  well  shown  in  Figure  3. 

The  Rayleigh  motion,  best  exhibited  on  the  vertical 
component,  is  a  wave  packet  with  almost  constant  frequency 
at  all  distances  at  which  it  has  been  observed.  This  does 
not  mean  that  the  medium,  known  to  be  layered,  is  not  dis¬ 
persive.  A  careful  correlation  of  peaks  and  troughs  shows 
that  dispersion  is  exhibited  as  a  progression  of  the  maximum 
amplitude  at  a  group  velocity  of  approximately  90  meters 
per  second  through  the  constant  frequency  wave  packet.  This 
is  demonstrated  in  Figure  5  in  which  the  vertical  component 
of  the  ground  motion  between  91  and  213  meters  is  plotted. 

The  transverse  motion  at  Suffield  has  been  identified  as  a 
Love  wave  on  the  basis  of  the  agreement  of  the  observed 
phase  velocity  curve  with  the  theoretical  curve.  (Figure  8). 
It  is  obvious  from  a  comparison  of  Figures  2  and  3  that  the 
long  dispersed  wave  train  of  Love  motion  was  not  observed 
along  the  profile  at  Suffield  as  it  was  at  the  Florissant 
test  site.  However,  it  has  been  established  (Kisslinger, 
e£  al#  1961;  Stauder  and  Kisslinger,  1962)  that  Love  waves 
are  generated  at  both  sites  with  an  asymmetrical  pattern  of 
radiation,  and  it  may  be  that  the  profile  at  Suffield 
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happened  to  lie  near  a  nodal  line.  The  Rayleigh  motion  at 
the  two  sites  differs  in  two  respects:  l)  the  motion  at 
Suffield  is  a  constant  frequency  wave  packet  whereas  the 
motion  at  Florissant,  though  still  exhibiting  a  fairly 
narrow  spectrum  (about  one  octave,  Figure  6),  is  a  normal¬ 
ly  dispersed  wave  train;  2)  the  vertical  component  dom¬ 
inates  the  Rayleigh  motion  at  Suffield  but  at  Florissant 
the  radial  component  is  bigger. 

A  conparison  of  the  Fourier  anqalitude  spectra  at  both 
sites  for  these  surface  waves  is  shown  in  Figure  6.  The 
Rayleigh  wave  spectrum  at  Florissant  is  peaked  at  approxi¬ 
mately  7  cps.  The  frequencies  at  which  the  amplitudes  are 
0.7  of  the  maximum  are  4.5  and  10  cps.  At  Suffield  the  ver¬ 
tical  spectrum  is  sharply  peaked  at  4  cps.  with  the  half 
power  points  at  3.5  and  5  cps.  Since  at  the  distances  in¬ 
volved  here  it  is  difficult  to  separate  the  beginning  of 
the  Rayleigh  wave  from  the  prograde  event,  each  event  affects 
the  spectrum  of  the  other.  This  is  most  noticeable  on  the 
radial  component  at  Suffield,  on  which  the  Rayleigh  wave  is 
relatively  small. 

The  theoretical  and  observed  dispersion  curves  for 
each  of  these  surface  waves  and  the  corresponding  , earth 
models  at  the  two  sites  are  shown  in  Figures  7  and  8. 

These  earth  models  were  selected  as  giving  the  best  agree¬ 
ment  with  the  observed  fundamental  Rayleigh  wave  curve 
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FIGURE  6.  FOURIER  SPECTRA,  FLORISSANT  AND  SUFFIELD 
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consistent  with  known  velocity  and  structural  data.  The 
curves  computed  from  these  models  for  the  other  types  of 
surface  waves  were  then  accepted  as  the  proper  theoretical 
curves. 

Since  the  method  of  analysis  yields  an  infinite  set  of 
possible  observed  phase  velocity  curves  (Sato,  1958,  Kiss- 
llnger,  i960),  it  is  necessary  to  select  one  of  these  on 
the  basis  of  independent  information.  At  both  sites,  P  and 
S  wave  velocities  have  been  directly  measured  and  layer 
thicknesses  are  known  from  refraction  surveys  and  logs  of 
shot  holes.  The  earth  models  used  for  the  theoretical 
curves  fit  well  with  the  known  velocities  and  layer  thick¬ 
nesses. 

An  interesting  feature  of  these  curves  is  that  the  Ray¬ 
leigh  wave  phase  velocity  is  greater  than  that  of  the  Love 
wave  over  part  of  the  frequency  band  of  interest.  An  in¬ 
spection  of  the  models  shows  that  the  shear  velocities  in 
the  upper  layers  at  both  sites  are  low.  In  the  clays  over- 
lying  the  more  rigid  substratum  the  value  of  Poisson* 3 
ratio  is  quite  high,  ranging  from  approximately  0.38  to 
0.49.  These  high  values  coupled  with  the  increase  of  the 
P-wave  velocity  at  the  first  boundary  without  a  correspond¬ 
ing  increase  in  the  shear  wave  velocity  cause  this  unusual 
relationship  of  Love  and  Rayleigh  phase  velocity  curves. 

The  high  values  of  Poisson’s  ratio  are  characteristic  of 
the  unconsolidated  surface  layers  and  the  increase  in  P-wave 
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velocity  without  similar  changes  in  shear  wave  velocity  is 
a  characteristic  of  shallow  water  tables. 

These  figures  also  show  the  dispersion  curves  of  the 
prograde  elliptic  events  discussed  previously.  At  Suffield 
the  vertical  and  horizontal  components  of  this  event  gave 
observed  phase  velocities  which  were  slightly  different. 

This  is  most  probably  due  to  the  interference  effects  with 
other  events  so  close  to  the  source.  However,  the  general 
agreement  of  these,  particularly  the  vertical  component, 
with  the  theoretical  curve  is  quite  good.  At  Florissant 
this  is  not  the  case.  As  already  mentioned  the  prograde 
motion  is  thought  to  be  a  combination  of  7  and  the  M21  mode, 
and  the  effect  of  this  interference  in  the  observed  curve 
is  not  known. 

Attenuation  coefficients,  oc  ,  were  determined  for  the 
individual  Fourier  components  of  the  fundamental  Rayleigh 
mode  from  the  Fourier  spectral  amplitude  at  varying  distances. 
The  values  were  determined  from  the  relation 

(5)  Ar  -  ^r-  e*P  (-  *  r)* 

where  Ar  -  Fourier  amplitude  at  distance  r, 

and  A0  •»  initial  amplitude. 

The  values  determined  for  the  two  test  sites  are  shown  in 
Table  2.  Because  the  variation  is  small  over  the  predomi¬ 
nant  frequency  range  of  the  Rayleigh  wave  at  both  sites, 
an  average  value  independent  of  frequency  was  assigned  to 
each  site  for  use  in  the  synthesis  calculations. 
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ATTENUATION  COEFFICIENTS 
Suffleld _ Florissant 


(°P8) 

(meters*-1) 

mm&nwmi 

■1  111  1  11 

3.5 

1.0037 

4.0 

0.0062 

4.0 

0.0044 

5.0 

0.0075 

4.5 

0.0040 

6.0 

0.0075 

5.0 

0.0046 

7.0 

0.0078 

5.5 

0.0057 

8.0 

0.0072 

9.0 

0.0078 

10.0 

0.0078 

11.0 

0.00B8 

Synthesis  of  Fundamental  Rayleigh  Mode 

The  synthesis  studies  from  these  data  are  In  two  categor¬ 
ies:  l)  the  attempt  to  synthesize  the  motion  in  the  source 
region,  i.e.  calculation  of  the  source  function,  and  2) the 
synthesis  of  motion  at  distances  sufficiently  far  from 
the  source  that  the  fundamental  Rayleigh  wave  has  separated 
from  the  other  events  such  as  direct  and  refracted  body 
waves  and  higher  modes. 

Synthesis  at  Distance.  The  motion  synthesized  at  135  and 
180  meters  frc  the  observed  motion  at  75  meters  at  Flor¬ 
issant  (Figure  2)  is  compared  with  the  observed  motion  in 
Figure  9.  Except  for  amplitude,  the  calculated  motion 
agrees  quite  well  with  the  observed  motion.  Slight  differ¬ 
ences  in  time  for  peaks  and  troughs  are  reasonably  accounted 
for  by  slight  inaccuracies  in  the  phase  velocity  data.  The 
smaller  amplitudes  on  the  synthesized  traces  are  largely 
due  to  the  effect  of  the  charge  size.  The  observed  motion 
at  75  meters,  the  input  to  the  synthesis  calculation,  was 


re¬ 


recorded  from  a  smaller  shot  than  that  which  produced  the 
observed  motion  at  135  and  180  meters*  Motion  recorded  at 
the  same  distance  from  both  of  these  shots  Indicated  that 
the  vertical  component  of  the  Rayleigh  motion  was  approx** 
lmately  1.3  times  greater  from  the  larger  shot.  The  mean 
value  of  the  ratio  of  observed  to  synthesized  at  135  meters 
Is  1.46  on  the  vertical  and  1.4  for  the  radial.  Figure  10 
Illustrates  the  agreement  of  amplitude  one  might  expect  when 
the  synthesized  motion  Is  calculated  from  motion  recorded 
from  the  same  shot.  This  figure  shows  the  vertical  motion 
synthesized  at  75  and  90  meters  from  observed  motion  at  90 
and  75  meters  respectively.  The  amplitude  agreement  Is 
quite  good.  Similar  calculations  for  the  radial  motion  pro¬ 
duced  the  same  results. 

The  synthesized  vertical  and  radial  motion  at  Suffield 
at  183  and  244  meters  as  calculated  from  the  observed  motion 
at  122  meters  is  shown  in  Figure  11.  Again  the  comparisons 
are  quite  favorable,  particularly  the  vertical  components. 
The  largest  motion  occurs  at  the  same  time  on  the  synthe¬ 
sized  and  observed  traces,  frequencies  are  the  same  and  the 
amplitudes  are  similar  and  die  off  rapidly  on  both  sides  of 
the  maximum  motion. 

The  observed  motion  from  183  meters  was  used  to  synthe¬ 
size  the  motion  at  122  meters  at  Suffield,  the  reverse  of 
the  process  discussed  above.  The  synthesized  motion  at  122 
meters  compared  very  favorably  with  the  observed.  This 
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sante  test  was  performed  at  Florissant  for  various  distance 
combinations  with  equally  good  results. 

From  these  examples.  It  Is  obvious  that  Sato's  tech¬ 
nique  of  synthesis  is  very  successful  in  predicting  or  re¬ 
producing  the  observed  motion  at  any  distance  from  the 
source  at  which  the  Rayleigh  wave  has  separated  from  the 
other  events  on  the  record.  This  implies  that  the  phase 
velocity  curves  must  be  reasonably  accurate,  that  attenua¬ 
tion  in  this  region  is  adequately  represented  by  the  assump¬ 
tions  made,  and  the  numerical  methods  are  satisfactory. 

Synthesis  in  the  Source  Region.  Before  examining  the 
calculated  source  functions  it  is  instructive  to  consider 
the  results  of  the  synthesis  at  Suffield  at  18  and  55  meters 
from  the  source.  The  motion,  shown  in  Figure  12,  was  cal¬ 
culated  from  the  observed  motion  at  122  meters. 

At  55  meters  the  maximum  calculated  vertical  motion  co¬ 
incides  with  the  maximum  observed  motion.  The  onset  of  the 
fundamental  Rayleigh,  especially  in  the  radial  component,  is 
somewhat  obscured  on  the  observed  trace  by  the  prograde 
event  preceding  it.  The  long  duration  of  the  oscillatory 
motion  observed  at  this  distance,  slightly  less  than  one 
wave  length  from  the  source,  is  especially  to  be  noted. 

At  18  meters,  approximately  one  quarter  wave  length 
from  the  source,  motion  is  observed  for  longer  than  one 
second.  The  frequency  of  this  extended  train  is  practically 


BS. 


0  1.0  2.0  0  1.0  2.0 

S  EC0N0S 

FIGURE  12.  OBSERVED  AND  SYNTHESIZED  MOTION,  SHORT  RANGE,  SUFFIELD 


-51- 


constant  and  crests  can  be  correlated  with  larger  amplitude 
phases  at  greater  distances  (Figure  5).  This  motion  is  re¬ 
produced  on  the  synthesized  traces.  A  comparison  of  the 
synthesized  and  observed  traces  reveals  that  the  maximum 
synthesized  vertical  motion  is  approximately  0.1  second 
later  than  the  maximum  motion  on  the  observed  trace.  The 
particle  motion  in  the  radial-vertical  plane.  Figure  13» in¬ 
dicates  that  the  fundamental  Rayleigh  motion  begins  at 
approximately  0.220  seconds,  which  corresponds  to  a  time 
0.040  seconds  before  the  observed  maximum  on  the  vertical 
trace.  However,  in  Figure  5  it  was  shown  that  the  maximum 
energy  is  propagating  at  a  group  velocity  of  90  meters  per 
second.  This  means  that  the  maximum  Rayleigh  motion  should 
occur  at  about  0.2  seconds.  The  interference  with  body 
waves  and  the  prograde  event  preceding  the  Rayleigh  wave 
alters  the  observed  traces  such  that  the  Rayleigh  onset  pre¬ 
dicted  from  the  particle  motion  diagram  would  not  be  accur¬ 
ate  but  would  give  the  time  at  which  the  Rayleigh  mode 
dominates  the  motion.  The  maximum  synthesized  motion  on 
the  radial  component  is  0.040  seconds  earlier  than  the 
vertical  maximum.  The  particle  motion  is  elliptic  retrograde. 
The  important  points  to  be  noted  here  are  the  delays  in  the 
calculated  maxima  and  the  long  duration  of  the  synthesized 
and  observed  traces,  which  are  particularly  well  shown  on 
the  vertical. 

The  motion  calculated  at  the  source  at  Suffield  from 
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the  observed  motion  at  122  meters  is  shown  in  Figure  14. 

The  most  striking  feature  is  the  long  oscillatory  nature 
of  the  source  function.  Figure  15  shows  a  comparison  of 
the  source  functions  calculated  from  different  distances 
and  different  shot  depths.  The  upper  two  traces  are  the 
vertical  component  calculated  from  a  shot  at  a  depth  of  5*5 
meters  and  recording  distances  of  122  and  183  meters  re¬ 
spectively.  The  lower  trace  was  calculated  from  motion 
observed  at  j.53  meters  from  a  shot  at  a  depth  of  2.3 
meters.  The  similarity  of  the  source  functions  is  obvious. 
This  shows  the  reproducibility  to  be  expected  from  calcula¬ 
tions  from  various  distances*  and  the  reproducibility  of 
the  motion  from  different  source  depths  in  the  same  medium 
for  these  cases  in  which  the  recording  distances  are  con¬ 
siderably  larger  than  the  source  depths. 

In  all  of  these  calculations  the  maxima  on  the  vertical 
and  radial  traces  are  approximately  0.3  seconds  after  origin 
time.  The  possibility  that  the  phase  velocity  curves  are 
grossly  in  error  is  ruled  out  when  one  considers  the  agree¬ 
ment  between  the  observed  and  calculated  motion  at  large 
distances*  aB  previously  shown  and  the  agreement  at  55 
meters.  Some  other  explanation  for  this  delay  must  be 
sought.  Before  pursuing  this  further,  let  us  consider  the 
source  functions  at  Florissant. 

Figure  16  shows  the  radial  and  vertical  source  functions 
which  were  calculated  from  the  observed  motion  at  75  and 
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FIGURE  15.  SOURCE  MOTION  FROM  DIFFERENT 
DEPTHS  AND  DISTANCES,  SUFFIELD 
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180  meters  from  the  source.  The  agreement  of  the  wave  forms 
on  the  radial  traces  is  excellent;  however,  there  is  a  dif¬ 
ference  in  the  wave  form  of  the  vertical  traces.  The  rela¬ 
tive  amplitudes  of  the  two  peaks  on  the  upward  aide  of  the 
zero  line  are  different.  The  observed  vertical  Rayleigh 
motion  at  75  meters.  Figure  2,  i3  3till  distorted  by  the  pro- 
grade  event  which  precedes  it.  The  accurate  determination 
of  the  Fourier  phase  angles  is,  therefore,  difficult,  and 
distortion  of  the  source  function  results.  These  events 
separate  as  the  distance  from  the  source  increases;  conse¬ 
quently,  the  source  functions  calculated  from  180  meters 
should  be  considered  more  reliable. 

The  source  functions  at  Florissant  are  oscillatory  in 
nature  but  much  shorter  in  duration  than  those  at  Suf field. 
The  maximum  amplitudes  occur  between  0.175  and  0.220  sec¬ 
onds  after  origin  time.  Here  also  the  comparisons  shown  at 
distance  (Figures  9  and  10)  give  one  confidence  that  the 
phase  velocity  data  used  in  the  calculations  are  reliable. 
This  means  that  we  must  look  to  another  source  for  the  de¬ 
lays  in  the  observed  maxima. 

In  the  series  of  tests  at  Suffield  a  500  pound  hemis¬ 
pherical  charge  was  detonated  on  the  free  surface  and  record¬ 
ed  between  153  and  366  meters  along  the  same  profile  as  the 
contained  shots.  The  Bource  functions  calculated  from  the 
observed  motion  at  213  meters  are  shown  in  Figure  17.  The 
calculated  motion  of  both  the  radial  and  vertical  components 
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ia  long  and  oscillatory  Just  as  the  motion  from  the  con¬ 
tained  shots.  However,  unlike  the  contained  shots,  the  max- 

irnum  amplitudes  occur  near  zero  time. 

Effects  of  Asymptotic  Approximation  In  Source  Region^ 

Cylindrical  spreading  was  assumed  In  the  calculation  of  the 


synthesized  motion.  The  pertinent  solution  of  the  wave 
equation  In  this  case  Is  the  Hankel  function  of  the  second 
kind.  For  distances  that  are  sufficiently  large  the  asymp¬ 
totic  expansion  of  the  Hankel  function  Is  used.  This  gives 
rise  to  the  r'i  factor  for  geometric  spreading.  It  was 
thought  that  the  use  of  the  asymptotic  approximation  In  the 
source  region  might  be  affecting  the  calculated  results.  To 
check  this  the  exact  value  of  the  amplitude  and  phase  of 
the  Hankel  function  for  each  frequency  was  calculated  for  dls- 
tances  of  18  meters  and  30  centimeters.  The  difference  be- 
these  computed  values  and  the  values  from  the  asymp¬ 


totic  expansion  were  used  to  determine  a  phase  and  amplitude 
correction  for  each  frequency.  The  radial  motion  at  Suffleld 
was  synthesized  using  the  corrected  values.  Comparisons  of 
the  synthesized  motions  at  18  meters  and  at  the  source  re¬ 
vealed  the  following!  1)  at  18  meters  there  Is  no  signifi¬ 
cant  change  In  the  motion:  amplitudes  and  frequencies  are 
the  same,  and  2)  at  the  source  there  Is  a  reduction  of 
amplitude  but  the  duration  of  motion  and  the  frequencies  are 
the  same.  Thus,  the  use  of  the  asymptotic  expansion  In  the 
source  region  seems  to  be  a  satisfactory  approximation  within 
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the  limits  of  precision  of  the  data  and  numerical  methods 
employed . 

Interpretation  of  Source  Functions,  In  this  discus¬ 
sion,  the  "source  function"  Is  the  motion  calculated  at  the 
surface,  over  the  shot  point,  which  would  produce  the  ob¬ 
served  Rayleigh  motion  at  distance  when  subjected  to  the 
given  dispersion.  For  a  near-surface  burst,  for  which  Ray¬ 
leigh  waves  are  produced  directly,  the  source  function  will 
take  the  form  of  the  elastic  response  of  the  medium  to  the 
applied  pressure  waveform.  For  a  shot  buried  so  deeply  that 
the  free  surface  does  not  participate  in  the  Initial  strong 
motion,  the  source  function  represents  the  response  of  the 
medium  to  a  point  force  applied  at  the  surface,  equivalent 
to  the  interaction  of  the  body  waves  from  the  explosion  with 
the  free  surface.  This  response  is  not  the  same  as  the 
displacement  of  the  boundary  of  the  equivalent  cavity  in  the 
conventional  model  of  an  explosive  source. 

The  accuracy  with  which  the  calculated  source  function 
represents  the  actual  behavior  at  the  epicenter  depends  on 
the  extent  to  which  the  input  spectrum  is  preserved  in  the 
observed  motion.  The  observed  spectrum  is  always  band- 
limited,  and  the  resulting  source  function  must  be  somewhat 
oscillatory.  A  sharply  peaked  spectrum  will  yield  a  highly 
oscillatory  source  function. 

Aki  (1960c)  interprets  the  source  function  on  the 
assumption  that  the  initial  phase,  the  same  for  all  fre- 
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quencies,  is  independent  of  the  layering  of  the  medium.  If 
the  distance  to  the  recorder  is  sufficiently  greater  than 
the  focal  depth,  as  it  is  for  all  of  the  data  used  here, 
then  the  initial  phase  of  the  source  function  depends  only 
upon  the  time  and  spatial  history  of  the  applied  force  at 
the  source.  Thus,  under  the  assumption  that  the  effect  of 
dispersion  due  to  the  layering  and  the  phase  changes  intro¬ 
duced  by  the  recording  instrument  are  accurately  removed  in 
the  process  of  phase  equalization,  the  polarity  and  symmetry 
of  the  source  function,  which  are  determined  by  the  initial 
phase,  can  be  interpreted  in  terms  of  idealized  forces  act¬ 
ing  at  the  origin. 

In  his  paper  Aki  compiled  idealized  source  functions 
for  the  vertical  component  associated  with  an  explosive 
step  and  an  explosive  impulse.  The  waveforms  (initial  phase) 
associated  with  these  sources  are  not  unique  to  them.  For 
example,  the  response  to  an  explosive  step  from  a  point 
source  is  the  same  as  that  from  a  horizontal  force  directed 
toward  the  receiver.  Figure  18  shows  the  source  functions, 
vertical  and  radial  components,  from  an  explosive  step  and 
an  explosive  impulse  applied  at  a  point  in  a  half-space. 

The  source  functions  calculated  for  Florissant,  Figure 
16,  closely  resemble  those  found  by  Aki  for  Rayleigh  waves 
from  earthquakes.  An  examination  of  Aki's  tabulated  source 
functions  and  the  present  results  indicate  that  the  assign¬ 
ment  of  these  derived  waveforms  to  a  particular  class 
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(even,  odd  or  some  combination  of  these)  depends  on  very 
small  departures  from  symmetry,  and  small  differences  in 
the  amplitudes  of  displacements  adjacent  to  the  maximum 
motion.  A  purely  even  or  purely  odd  function  can  only  re¬ 
sult  from  a  two-dimensional  source.  Since  the  source  in 
these  experiments  is  known  to  be  a  point,  the  selection  of 
an  even  or  odd  function  was  excluded,  even  though  this  re¬ 
quires  placing  confidence  in  minor  variations  in  the  wave¬ 
form.  If  the  interpretation  is  constrained  to  a  three- 
dimensional  solution,  then  the  result,  based  on  the  record 
at  180  meters,  for  which  the  second  peak  on  the  vertical  is 
greater  than  the  first,  but  smaller  than  the  trough,  while 
the  trough  on  the  radial  component  is  slightly  smaller  than 
the  following  crest,  agrees  best  with  an  explosive  step. 

The  results  for  75  meters  indicates  a  down  pulse  in 
force  rather  than  an  explosive  3tep.  Since  it  is  known  that 
motion  other  than  the  fundamental  Rayleigh  wave  is  included 
in  the  record  at  this  station,  little  confidence  can  be 
placed  in  the  small  differences  between  the  two  results 
that  require  the  different  interpretations. 

It  must  be  recognized  that  there  is  no  reason  to  demand 
that  the  driving  force  in  the  field  experiments  be  either 
impulsive  or  a  step  function,  or  that  the  initial  phase 
of  all  Fourier  components  be  the  same.  The  effect  of  a 
driving  force  that  decays  monotonically  or  oscillates  would 
be  to  destroy  the  ideal  symmetry  and  extend  the  source 
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function  in  time. 

The  source  functions  from  the  contained  shots  at  Suf- 
field,  Figures  14  and  15,  are  long  in  duration  and  of  almost 
constant  frequency,  a  consequence  of  the  narrow  frequency 
spectrum  (Figure  6).  These  source  functions  cannot  be  in¬ 
terpreted  as  those  at  Florissant,  because,  as  Aki  states 
(1960c,  p.  2408):  "If  the  spectrum  is  peaked,  the  source 
function  will  take  an  oscillatory  form,  which  is  undesir¬ 
able  for  the  clear  identification  of  its  type." 

Before  considering  why  the  frequency  spectrum  is  so 
narrow  at  Suffield  let  us  examine  the  source  functions  from 
the  500  pound  surface  shot.  Figure  17.  As  in  the  case  of 
the  contained  shots,  the  waveforms  of  both  components  oscil¬ 
late  with  constant  frequency,  and  the  spectrum  of  the  verti¬ 
cal  component  peaks  at  4  cps  with  the  same  band  width.  The 
fact  that  the  peak  of  the  spectrum  does  not  change  with  a 
fifty-fold  increase  in  charge  size  requires  that  some  factor 
other  than  total  input  energy  is  controlling  the  spectrum. 
The  spectrum  for  a  given  medium  can  be  affected  by  crater  or 
cavity  size  and  by  the  layering  present.  Since  charge  size 
has  no  appreciable  effect  here,  the  near  surface  layering 
must  dominate  the  spectrum.  The  contained  shots  and  the 
surface  shot,  two  distinctly  different  sources,  produce 
motion  at  distance  which  is  identical  in  the  frequency  of 
the  waveform  but  differ  in  phase  by  TT  .  This  leads  to  the 
conclusion  that  the  layering  is  important  in  the  Rayleigh 
wave  generation  process  at  this  site. 
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It  was  shown  above  (Figure  12)  that  constant  frequency 
motion  was  observed  for  longer  than  one  second  at  a  dis¬ 
tance  approximately  one-quarter  wave  length  from  the  source. 
These  observations  suggest  that  energy  Is  being  transmitted 
from  the  source  region  for  a  period  of  time  much  longer 
than  the  time  duration  of  the  explosion  Itself.  This  Is 
typical  of  a  tuned  or  resonant  mechanical  oscillator  with 
low  damping  and  leads  to  the  belief  that  the  Rayleigh  motion 
at  this  site  is  produced  by  a  resonance  phenomenon. 

There  Is  additional  Information  to  support  this  view¬ 
point.  In  addition  to  the  contained  ten-pound  shots  fired 
at  Suffleld,  ten  pounds  was  fired  on  the  surface  and  3 


meters  above  the  free  surface.  Neither  of  these  shots  gen¬ 
erated  any  Rayleigh  waves.  In  August,  1961,  a  100-ton  sur¬ 
face  burst  was  fired  at  this  test  site.  In  reviewing  the 
seismograms  from  this  test,  Klsslinger  (1963)  pointed  out 
that  the  Rayleigh  wave  generated  by  smaller  shots  was  not 
present.  A  review  of  records  from  earlier  trials  at  this 
teat  site  showed  Rayleigh  waves  with  the  same  period  and 
group  velocity  as  those  from  the  smaller  charges  to  be  ab¬ 
sent  from  all  of  the  surface  shots  greater  than  five  tons. 
Retrograde  elliptic  motion  following  the  air  blast  arrival 
is  found  on  seismograms  from  these  large  shots  but  the 
arrival  times  are  much  too  early  for  this  to  be  the  same 

event  as  that  recorded  from  the  smaller  shots.  Klsslinger 
suggested  that  the  100-ton  shot  acted  as  a  force  downward 
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on  the  lower  layers  and  that  the  longer  period  Rayleigh 
motion  resulting  was  controlled  by  the  higher-velocity 
deeper  layers.  The  size  of  the  crater  for  the  100-ton 
shot  would  tend  to  support  this  conclusion.  The  entire 
thickness  of  the  low  velocity  surface  layer  was  involved 
in  the  source  region. 

The  evidence  presented  here  leads  one  to  the  conclu¬ 
sion  that  the  response  of  the  free  surface  to  the  explosive 
loading  that  results  in  the  production  of  the  fundamental 
Rayleigh  waves  at  the  two  sites  is  entirely  different.  At 
Florissant  the  source  function  appears  to  be  independent  of 
the  layering  whereas  at  Suffield  it  appears  to  be  dominated 
by  the  near-surface  layering.  For  small  charges  at  Suf¬ 
field  this  effect  does  not  depend  on  whether  the  charge  is 
at  the  surface  or  buried. 

It  was  noted  in  a  previous  section  that  there  is  a 
time  delay  for  the  maximum  amplitude  of  the  source  functions 
from  the  contained  explosions  which  did  not  appear  in  the 
source  functions  from  the  surface  shot.  This  suggests  that 
there  is  a  time  delay  in  the  generation  process  of  Rayleigh 
waves  from  a  contained  explosion  which  does  not  exist  for 
the  surface  source.  This  time  delay  is  longer  than  P  or  S 
wave  travel  times  from  the  source  to  the  surface.  This 
point  merits  further  study.  The  data  from  a  small  number 
of  tests  under  nearly  similar  conditions  are  not  sufficient 
to  allow  any  definite  conclusion. 
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Rmirce  Functions  from  Higher  Modes. 

The  presence  of  prograde  elliptic  motion  preceding  the 

fundamental  Rayleigh  mode  at  both  teat  sites  has  already 
been  discussed.  Source  functions  were  calculated  for  these 


events . 

Figure  19  shows  the  radial-  and  vertical-component 
source  functions  calculated  from  the  observed  motion  at  180 
meters  from  Florissant.  The  wave  forms  are  oscillatory  and 
more  complex  than  were  the  fundamental.  Source  functions 
calculated  from  various  distances  were  not  reproducible. 
This  suggests  that  the  phase  velocity  curve  used  for  the 


synthesis  of  the  Fourier  components  Is  not  correct.  This 
is  not  surprising  In  view  of  the  Interpretation  of  this  pro 
grade  motion  as  a  combination  of  7  and  the  mode  of  the 
Rayleigh  wave. 


An  interpretation  of  these  source  functions  In  terms  of 
a  force  system  cannot  be  justified  In  view  of  the  lack  of 
reproducibility  observed.  If  the  slope  of  the  phase  veloc¬ 
ity  curve  is  m  error  or  If  the  values  of  the  velocity  at 
each  frequency  are  In  error  by  some  value,  the  resulting 
calculations  will  be  in  error.  There  Is  good  reason  to  sus¬ 


pect  that  this  is  tho  case. 

The  motion  of  the  H21  branch  observed  at  183  meters  at 
Suffleld  was  used  to  calculate  the  source  functions  shown 
in  Figure  20.  This  motion  was  calculated  using  a  phase 
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velocity  curve  which  was  slightly  faster  than  the  observed 
vertical  curve  in  the  frequency  range  from  6  to  11  cps.  and 
slightly  lower  than  the  observed  radial  curve  in  the  same 
range.  This  curve  was  used  because  tests  with  the  observed 
curves  indicated  that  both  were  somewhat  in  error.  The  mo¬ 
tion  synthesized  at  the  source  using  the  observed  curves 
was  very  similar  to  that  using  the  average  curves  but  better 
time  agreement  between  the  radial  and  vertical  components 
was  attained  with  the  average  curve. 

The  source  function  is  much  simpler  than  that  from  the 
fundamental  Rayleigh  wave  at  this  situ.  The  spectrum  for 
this  event.  Figure  6,  is  broader  than  that  for  the  fundamen¬ 
tal  and  peaks  at  a  higher  frequency.  The  ground  motion  for 
this  event,  shown  in  Figure  3,  remains  pulse-like  ir.  nature 
over  the  range  of  recording  and  does  not  exhibit  strong  dis¬ 
persion. 

The  difference  in  the  nature  of  the  source  function  be¬ 
tween  the  M21  mode  and  the  fundamental  can  be  interpreted 
to  mean  that  the  generation  process  is  somewhat  different 
for  the  two  events.  The  Bource  function  for  the  M21  mode 
can  be  interpreted  UBing  Aki's  criterion.  The  waveform  of 
the  source  function  is  similar  to  that  associated  with  a 
down  pulse  propagating  in  three  dimensions.  This  is  some¬ 
what  surprising  in  view  of  the  fact  that  the  explosion  was 
contained.  This  event  could  not  be  identified  for  the  500 
pound  shot,  because  its  expected  time  of  arrival  coincided 
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v;ith  the  air  blast  arrival  and  the  large  motion  which  fol¬ 
lows  . 

The  simple  nature  of  this  source  function  strengthens 
the  view  that  something  other  than  the  shot  is  controlling 
the  generation  of  the  fundamental  Rayleigh  motion  at  Suf- 
field. 

Synthesis  of  Love  Waves, 

The  fact  that  explosions  produce  prominent  horizontal 
transverse,  SH,  motion  has  been  discussed  in  detail  for  the 
Florissant  test  site  (Kisslinger,  et  al,  19&1  and  Chapter  V). 
The  long  period  motion  observed  on  the  transverse  trace  has 
been  identified  as  a  Love  wave.  The  dispersion  curve  is 
shown  in  Figure  7  and  the  Fourier  amplitude  spectrum  in 
Figure  6. 

The  Love  motion  from  a  profile  at  Florissant  was  used 
to  calculate  the  source  functions  shown  in  Figure  21.  The 
reproducibility  of  the  wave  form  from  the  two  distances  from 
which  calculations  are  shown,  150  and  180  meters,  is  ex¬ 
cellent.  The  waveform  is  oscillatory  but  basically  fairly 
simple.  The  duration  and  shape  are  similar  to  those  from 
the  fundamental  Rayleigh  at  this  site.  The  frequency  of  the 
waveform  is  approximately  10  cps  which  corresponds  to  the 
peak  of  the  Fourier  spectrum.  Clockwise  motion  (down  on 
the  trace)  about  the  source  is  indicated  for  this  azimuth. 
Conclusions. 


1.  At  distances  at  which  the  fundamental  Rayleigh 
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motion  is  separated  from  the  faster  events,  the  synthesized 
Rayleigh  motion  matches  the  observed  motion  quite  closely. 
This  indicates  that  the  method  of  synthesis  of  Fourier  com¬ 
ponents  is  useful,  for  one  can  accurately  predict  the 
motion  at  any  distance  as  long  as  the  dispersive  properties 
do  not  significantly  change,  and  the  observed  motion  does 
not  represent  an  interference  of  two  or  more  events. 

2.  The  source  functions  of  the  fundamental  Rayleigh 
waves  from  the  two  test  sites  are  quite  different.  Both 
are  oscillatory  in  nature  but  the  very  narrow  spectrum  at 
Suffield  results  in  a  long  duration  constant  frequency 
waveform.  This  has  been  interpreted  to  mean  that  there  is 

a  distinctly  different  response  of  the  near-surface  materials 
to  the  explosion  at  the  two  sites.  The  evidence  at  Suffield 
leads  to  the  conclusion  that  a  resonance  phenomenon  con¬ 
trolled  by  the  near  surface  layering  strongly  affects  the 
character  of  the  Rayleigh  wave,  whereas  at  Florissant  the 
layering  plays  no  comparable  part  in  the  generation  process. 

3.  The  maximum  amplitudes  of  the  source  functions  from 
contained  explosions  are  delayed  as  long  as  0.3  seconds  after 
the  origin  time.  This  is  not  the  case  for  a  surface  shot 

at  Suffield.  The  conclusion  is  that  the  delay  is  due  to 
the  generating  mechanism  for  Rayleigh  waves  from  a  contained 
explosion,  requiring  that  energy  proceed  from  the  source  re¬ 
gion  for  a  relatively  long  time  after  the  detonation. 
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4.  Calculations  of  source  functions  for  Love  waves  at 
Florissant  from  different  distances  along  one  azimuth  were 
very  similar. 

5.  The  classification  of  calculated  source  functions 
into  particular  types  in  order  to  interpret  the  character 
of  the  input  force  is  difficult.  Confidence  must  be  placed 
in  minor  differences  in  the  calculated  waveforms,  differences 
that  may  well  be  within  the  acouracy  of  the  data  and  the 
analytic  procedures. 


CHAPTER  IV 


RADIATION  PATTERNS  PROM  RIPPLE-FIRED  SHOTS 

Background  and  References. 

For  many  years  quarry  and  mine  operators  have  been 
utilizing  ripple -fired  shots  to  minimize  the  damaging  ef¬ 
fects  of  their  blasting.  In  practice,  the  ground  amplitude 
or  particle  velocity  are  reduced  by  a  trial  and  error  de¬ 
termination  of  the  optimum  spacing  and  delay  times.  Willis 
(1963)  studied  velocity  spectra  from  instantaneous  and  .017 
sec.  delayed  multi-hoie  shots.  He  found  amplitude  attenua¬ 
tion  from  ripple-firing  to  be  considerably  more  for  body 
waves  than  surface  waves.  The  effect  is  dependent  on  fre¬ 
quency  spectrum  and  wave  length  of  the  seismic  signal. 

Pollack  (1963)  observed  the  predicted  relationship  be¬ 
tween  the  spectrum  of  a  single  shot  and  that  of  the  same 
signal  summed  with  time  delays  in  a  series  of  controlled 
ripple-fired  quarry  blasts.  He  assumed  a  smooth,  damped 
oscillatory  waveform  of  approximately  .001  sec.  rise  time 
and  .004  sec,  duration  (simple  body  wave  arrival)  as  the 
unit  input  signal  and  compared  observed  and  predicted  spec¬ 
tra  for  3,  7  and  15  holes  delayed.  In  total,  .00166, 

.00375,  and  .00664  sec.  The  effect  of  any  time  delays  due 
to  spatial  distribution  on  the  shots  was  not  considered. 
These  would  introduce  an  azimuthal  variation  into  the  re¬ 
sultant  signal. 

Ripple-firing  effects  on  surface  waves  are  complicated 
by  the  addition  of  a  frequency-dependent  velocity,  so  that 
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the  equivalent  time  delay  corresponding  to  a  spatial  separa¬ 
tion  depends  also  on  frequency. 

As  a  part  of  an  overall  study  of  controlled  radiation 
patterns  of  P,  SH,  and  Rayleigh  waves  from  explosions  a 
series  of  field  experiments  in  ripple-firing  was  conducted 
in  the  summer  of  1962  (shot  points  65-72,  Table  1).  This 
section  of  the  report  presents  the  results  of  the  experiments. 
Data. 

Several  types  of  shot  point  patterns  and  various  time 
delays  were  used  in  the  experiments  in  order  to  test  the 
effects  of  the  spatial  distribution  of  the  individual  charges 
making  up  the  shots,  as  well  as  the  eifects  of  varying  the 
delay  times.  Motion  was  recorded  30  meters  from  the  center 
of  the  shot  pattern  at  s oc  equally  spaced  azimuths  with 
Sprengnether  three-component  seismographs  and  the  12  SIE 
vertical  geophones  were  spaced  between  the  Sprengnethers . 

The  basic  separation  unit  between  charges  in  the  shot  pat¬ 
terns  was  1.5  meters  and  the  unit  charge  size  was  0.25 
pound  of  dynamite'.  In  the  subsequent  discussion  time  de¬ 
lays  are  given  in  the  form  of  the  effective  velocity  of 
propagation  of  the  detonation. 

In  addition  to  calibration  shots  (S.P.  69a,  b)  of  1.25 
and  0.25  peund  at  the  center  of  the  detector  array,  ten 
multi-hole  shots  were  detonated  in  the  series* 

3. P.66?  5  shots  simultaneously  in  a  vertical  hole 

7.6  meters  deep 
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S.P.  65: 

S.P.  67a: 
S.P.  67b: 

S.P.  68a: 

S.P.  68b: 

S.P.  70: 


5  shots  detonated  upward  at  200  raeters/sec. 
in  a  vertical  hole  7*6  meters  deep. 

5  shots  simultaneously  in  a  horizontal  line 

5  shots  at  210  meters/sec.  in  a  horizontal 
line 

5  shots  at  103  meters/sec.  in  a  horizontal 
line 

5  shots  at  82  meters/sec.  in  a  horizontal 
line 

2  shots  at  190  meters/sec.  in  a  horizontal 
line 


S.P.  71:  2  shots  at  91  meters/sec.  in  a  horizontal 

line 

S.P.  72a:  8  shots  in  two  horizontal  lines  (couple)  of 

4  charges  detonated  in  opposite  directions 
at  185  raeters/sec. 

S  P  72b:  8  shots  in  two  horizontal  lineB  (couple)  of 

4  charges  detonated  in  opposite  directions 
at  91  meters/sec. 

All  delayed  shot3  were  fired  with  a  sequential  timer  utiliz¬ 
ing  rotating,  staggered  cams,  capable  of  firing  up  to  15 
shots.  Firing  times  of  the  first  and  last  shots  were  re¬ 


corded  on  the  SIE  timing  trace  for  determining  detonation 


velocities. 

Figures  22  through  28  present  the  reduced  seismograms, 
showing  the  radiation  patterns  produced  by  the  various  hori¬ 
zontal  shot  patterns.  S.P.  68b  is  omitted  because  a  malfunc¬ 
tion  of  the  timer  resulted  in  an  unusually  large  time  lag  in 
the  s.ot  sequence.  AI30,  in  the  case  of  S.P.  67a#  it  is 
doubtful  that  the  total  charge  detonated.  Azimuths  are  all 
measured  clockwise  from  the  northward  end  of  the  pattern. 
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GROUND  MOTION,  5  HOLE  PATTERN 


FIGURE  24.  GROUND  MOTION,  5  HOLE  PATTERN 
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URE  25.  GROUND  MOTION,  2  HOLE  PATTERN 


180 


FIGURE  26.  GROUND  MOTION,  2  HOLE  PATTERN 


FIGURE  27,  GROUND  MOTI  ON,  COUPLE 
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Data  Analysis. 

It  is  evident  from  the  figures  that  the  temporal  and 
spatial  distributions  of  these  shots  have  severe  effects  on 
the  radiation  patterns.  It  is  also  apparent,  from  the  sim¬ 
ilarity  of  records  at  points  symmetric  with  respect  to  the 
propagating  source,  that  there  are  no  significant  anomalies 
introduced  by  differences  in  path  characteristics. 

The  analysis  procedure  adopted  for  these  data  is  an 
attempt  to  explain  the  radiation  patterns  by  simple  superpo¬ 
sition  of  the  unit  seismic  signal  through  Fourier  synthesis. 

If  (t)  is  the  recorded  trace  from  the  unit  charge 
from  which  the  delayed  system  will  be  constructed,  then 

is  the  signal  from  N  shots  delayed  after  the  first  by  time 
A  tn,  where  A  tn  is  the  combined  delay  due  to  both  temporal 
and  spatial  separations,  i.e. 

=  Dn  +  -A  Rn 

where: 

Dp  =  firing  time  delay, 

ARn  «  additional  path  from  nth 
charge  to  receiver 

V  =  velocity  of  propagation  of  signal, 
possibly  frequency  dependent. 

In  terms  of  spectra,  if 

P1  <")  °  l  r,  (t)e-lwtdt 

-m 
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is  the  "unit"  spectrum,  then  for  the  composite  signal 

Jj  Za  e'lw tdt 
-  JjiWlutat 

„  e~iujAtn  Px  (to), 

If  the  amplitudes,  but  not  waveforms  of  the  individual 
shots  vary,  the  summation  contains  an  amplitude  factor.  ThiB 
is  the  case  for  the  vertical  line  of  charges,  where  such  a 
factor  was  used  in  the  computation.  Then  the  composite 
spectrum  is  given  by 

F(6J)  =  n  e“luj4tn  F1  (<w) 
where  the  individual  charges  give  signals 

fn  <*>  =  Vl  <*-***>• 

Summarizing,  the  spectrum  of  the  composite  signal  is 
given  by  the  spectrum  of  the  unit  signal  modified  by  the 
multiplicative  function: 

EA  e-i^£tn  for  variable  amplitudes  and 
n  delays, 

y^e"1  ^  Atn  for  variable  delays  only 
y~y.lcunat 


for  constant  delays,  A  t,  and  no 
spatial  effect. 


For  fixed  firing-time  delays,  the  variation  in  the  total 
delays  is  introduced  by  the  spatial  separation  of  the  indi¬ 
vidual  charges.  Without  this  factor  there  is  no  azimuthal 
effect  inherent  in  time  delays. 
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In  the  analysis,  synthetic  seismograms  were  computed 
by  performing  the  inverse  Fourier  transform  operation  on  the 
computed  F  (oj).  In  the  most  general  case  of  variable  de¬ 
lays  and  charges  the  synthetic  seismogram  is 

f3(o  =  1/  /[[v'1-tn]  *lwt  Fi (lu)  d“'- 

Figure  29  shows  the  nature  of  the  unit  signal  used  in 
this  procedure.  The  vertical  and  radial  components  and  the 
moduli  of  the  Fourier  Transforms  are  presented.  ThiB  motion 
was  produced  at  30  meters  by  a  calibration  shot  of  0.23 
pound  of  dynamite.  Both  the  vertical  and  radial  components 
of  motion  were  found  to  be  very  similar  in  all  azimuths  and 
one  record  was  arbitrarily  chosen  as  the  standard,  f^  (t). 
Subsequent  figures  show  theoretical  seismograms  for  several 
of  the  experimental  shot  patterns . 

Because  the  predominant  energy  in  the  standard  signal 
appears  to  be  the  fundamental  mode  Rayleigh  wave,  the  cor¬ 
responding  frequency-phase  velocity  relationship  was  used  in 
calculating  the  equivalent  time  delays  due  to  spatial  dis¬ 
tributions.  The  dispersion  curve  is  well-known  for  the  ex¬ 
perimental  area,  (see  Chapter  III). 

Figure  30  shows  the  theoretical  seismograms  for  the 
vertically  stacked  shots.  In  this  case  there  is  no  azimuth¬ 
al  variation  of  the  signal.  Because  of  the  amplitude-depth 
function  used,  the  upper  shot  controls  the  appearance  of  the 
seismogram  resulting  in  little  apparent  change  of  the  wave¬ 
form.  This  depth-dependence  function  was  constructed  from 
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observation  of  amplitude  change  over  the  same  depth  range 
for  other  shot  points  in  the  field  program.  Table  3>  in 
which  unit  amplitude  is  assigned  to  the  shallowest  shot,pre- 
sents  this  relationship.  Zero  time  corresponds,  in  the  case 
of  the  delayed  shot,  to  the  firing  time  of  the  top  charge. 

The  actual  seismograms  are  those  from  180°  azimuth.  The 
high  frequency  (50  to  80  cps)  body  wave  events  at  the  be¬ 
ginning  of  the  record  are  not  present  on  the  theoretical 
seismograms  because  of  the  50  cps  cut-off  on  the  standard 
spectrum,  a  result  of  digitalizing  increment .  In  addition, 
the  anplitude-depth  function  used  does  not  apply  to  the  body 
waves,  having  been  determined  for  the  Rayleigh  wave. 

Table  3 

AMPLITUDE  -  DEPTH  RELATIONSHIP  FOR  VERTICAL  PATTERNS 

A 

1.00 
.30 
.15 
.10 

.08 

Figures  31  and  32  show  theoretical  seismograms  for  the 
five-shot  horizontal  line  fired  simultaneously  and  at  108 
meters/sec . 

Figure  33  ehowB  theoretical  seismograms  for  the  two-shot 
horizontal  line  fired  at  91  meters/sec.  With  the  exception 
of  the  high  frequency  components  in  the  records,  the  agree¬ 
ment  is  generally  good.  On  the  vertical  components,  traces 
were  mistakenly  synthesized  at  azimuths  midway  between  re- 
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FIGURE  31.  THEORETICAL  GROUND  'MOTION,  5  HOLE  PATTERN 
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cording  azimuths.  However,  the  change  is  gradual  enough 
that  the  comparison  can  be  made.  The  shot  at  190  me ters/sec. 
is  similar  except  that  higher  speed  reduces  the  severity  of 
azimuthal  variations. 

Figure  34  is  the  theoretical  pattern  for  the  double 
horizontal  line  at  91  meters/sec.  No  records  were  obtained 
at  240°  and  300°,  which  should  be  equivalent  to  60°  and 
120°  by  symmetry.  The  185  meters/sec.  pattern  is  very  sim¬ 
ilar.  In  this  case  agreement  is  generally  good,  although 
there  are  some  significant  amplitude  differences. 

Summary  and  Conclusions. 

While  there  are  discrepancies,  probably  due  to  loss  of 
high  frequency  components,  incomplete  charge  detonations, 
and  errors  in  assignment  of  velocities  to  frequency  compo¬ 
nents,  radiation  patterns  observed  from  ripple-fired  shot 
arrays  can  be  reasonably  well  explained  theoretically  by 
simple  superposition.  Fourier  theory  indicates  that  a  per¬ 
turbation  to  the  unit  spectrum  in  the  form  of  the  sum  of  a 
product  of  amplitude  and  phase  change  factors  yields  the 
composite  spectrum.  The  total  time  delay  consists  of  a 
fixed  increment  and  a  spatial  time  equivalent,  depending 
upon  receiver  position.  The  spatial  distribution  gives  rise 
to  the  azimuthal  radiation  pattern.  Synthetic  seismograms 
were  computed  and  plotted  on  an  IBM  1620  for  all  recording 
azimuths.  Comparison  with  observed  motion  indicates  gen¬ 
erally  good  agreement. 
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The  severity  of  the  azimuthal  variations  that  can  be 
produced  suggest  that  maxima  or  minima  of  Rayleigh  wave 
energy  from  explosions  can  easily  be  oriented  if  ripple¬ 
firing  is  feasible,  and  if  the  shot  point  layout  can  be 
arbitrarily  selected. 
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CHAPTER  V 

i 

HORIZONTALLY  POLARIZED  SHEAR  WAVES 


Introduction. 

The  properties  of  shear  waves  generated  by  underground 
explosions  was  the  first  subject  of  investigation  in  this  re¬ 
search  program.  Simple  methematlcal  models  of  an  explosion 
as  a  source  of  seismic  waves  fail  to  predict  the  existence 
of  shear  waves  proceeding  directly  from  the  neighborhood  of 
the  explosion.  The  observed  presence  of  shear  waves  polar¬ 
ized  in  the  plane  of  incidence  (SV  waves)  is  readily  ex¬ 
plained  by  the  conversion  of  congressional  energy  at  subsur¬ 
face  interfaces  and  at  the  free  surface  of  the  earth,  and  so 
their  occurrence  cannot  be  taken  as  a  significant  departure 
from  simple  explosion  theory.  Horizontally  polarized  shear 
waves  (SH  waves)  cannot  be  generated  in  this  manner,  and  so 
it  was  decided  that  the  best  way  in  which  to  study  explosion¬ 
generated  shear  waves  was  to  concentrate  on  SH  waves.* 

The  existence  of  SH  motion  of  large  amplitude  on  seis¬ 
mograms  from  explosions  had  long  been  known  to  seismologists 
working  with  three -component  records  of  quarry  blasts,  but 
these  had  never  been  studied  systematically.  It  was  deter¬ 
mined  very  quickly  from  the  first  experiments  carried  out  in 
this  program  that  SH  motion  is  always  produced  by  buried 

*It  is  possible  for  apparent  SH  motion  to  arise  by  reflection 
and  head  wave  propagation  from  subsurface  beds  that  are  not 
parallel  to  the  earth1 s  surface. 
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shots,  and  that  these  waves  have  some  well-defined  proper¬ 
ties.  A  paper,  hereafter  referred  to  as  Paper  I,  in  which 
these  properties  are  explained  and  possible  generating  mech¬ 
anisms  are  discussed  was  published  in  1961  (Kisslinger,  £t 
al,  1961).  The  important  properties  may  be  summarized  as 

follows. 

1)  Both  body  shear  waves  and  Love  waves  are  generated  by  the 
explosion.  These  waves  have  frequencies  in  the  same  band 
as  the  corapressional  and  Rayleigh-type  waves. 

2)  Love  waves  are  propagated  with  a  pronounced  azimuthal  var¬ 
iation  in  amplitude,  marked  by  a  pattern  of  nodal  lines 
with  reversed  polarity  on  either  side  of  the  nodes. 

3)  The  observed  reversals  in  polarity  are  not  apparent  re¬ 
versals  of  the  phase  associated  with  the  maximum  ampli¬ 
tude  as  a  result  of  non-symmetry  of  the  dispersive  proper¬ 
ties  of  the  layered  system.  Reversal  of  the  entire 
horizontal  transverse  waveform,  down  to  fine  details,  has 
been  observed  in  many  cases,  and  the  pattern  of  radiation 
has  been  found  to  be  maintained  over  a  range  of  distances. 

4)  In  soft  materials,  such  as  clay  or  alluvium,  the  body 
shear  waves  are  roughly  one-fourth  to  one-third  the  ampli¬ 
tude  of  the  corapressional  waves,  and  Love  waves  have 
average  amplitudes  about  the  same  fraction  of  the  Ray¬ 
leigh  wave,  although  they  may  be  almost  as  big  at  the 
maximum  of  a  lobe.  In  these  materials  the  azimuthal 
distribution  of  radial  and  vertical  motion  is  completely 
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symmetrical,  and  the  SH  motion  is  undoubtedly  independent 
of  these. 

5)  In  brittle  materials,  such  as  sandstone  and  limestone, 
all  types  of  motion  are  somewhat  asymmetrical,  and  the 
shear  waves  may  be  as  large  as  the  other  motion  recorded. 
A  nuclear  shot  in  granite  has  confirmed  the  large  ampli¬ 
tudes  of  shear  waves  from  shots  in  competent  materials. 

6)  The  shear  energy  that  is  directly  produced  by  the  shot 
almost  certainly  is  formed  in  the  transition  region  of 
non-linear  behavior  around  the  explosion.  Processes  in 
this  region  are  only  partly  understood,  but  recent  work 
has  made  significant  progress  toward  clarifying  them 
(Bishop,  1962).  The  present  investigators  prefer  tne 
growth  of  radial  cracks  as  an  adequate  source  of  shear 
waves  with  the  observed  properties. 

Basic  Observations. 

The  decision  was  made  at  the  beginning  of  the  study  to 
put  aside  all  preconceived  notions  of  what  properties  ex- 
plosion-generated  shear  waves  should  have,  and  of  how  they 
might  be  generated.  Experiments  were  carried  out  that  were 
designed  to  provide  basic  data  on  these  questions.  Since 
the  first  experiments  were  completed,  hundreds  of  additional 
seismograms  have  been  collected  in  connection  with  other 
aspects  of  the  research  that  contain  further  evidence.  All 
of  the  information  supports  the  original  conclusions  as 
published  in  Paper  I,  and,  indeed,  strengthens  the  arguments 
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present  ed  there . 

The  fundamental  point  la  that  for  a  point  aource  In  a 
site  composed  of  essentially  horizontal  beds  of  uniform  com- 
position,  the  amplitude  of  the  SH  motion  la  different  In 
different  directions,  and  most  Important  of  all,  the  polar¬ 
ity  of  the  waveform  Is  completely  reversed  In  some  direc¬ 
tions  compared  to  others.  The  Importance  of  this  fact 
lies  in  the  difficulty  that  would  have  been  encountered  In 
trying  to  explain  the  origin  of  this  motion  If  It  had  been 
found  to  be  radiated  uniformly  In  all  directions.  Since 
numerous  observations  confirm  that  the  radial  and  vertical 
components  of  motion  are  uniform  with  azimuth  under  the 
same  conditions,  there  can  be  no  doubt  that  the  mechanism 
of  veneration  of  the  close-ln  SH  motion  Is  quite  independent 


of  the  other  modes  of  propagation. 

A  convenient  manner  In  which  to  portray  the  amplitude  data 
IS  by  means  of  polar  plots  of  the  amplitude  of  a  selected 
phase,  one  which  can  be  correlated  at  all  Instruments.  If 
a  similar  analysis  were  attempted  using  dispersed  Love 
waves  recorded  at  different  distances  and  over  different 
paths,  It  would  be  necessary  first  to  phase  equalize  the 
signals,  using  appropriate  phase  velocity  curves  and  one  of 
the  techniques  of  Chapter  III,  before  attempting  to  deter¬ 
mine  the  proper  polarity  to  assign  to  each  record.  However, 
in  these  records.  In  which  the  layering  and  distances  are 
Identical,  It  Is  possible  to  pick  a  prominent  phase  and 
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correlate  It  around  the  circle. 

Examples  of  polar  plots  for  Florissant  are  given  In 
Paper  I,  and  data  for  Suffield,  Alpha  and  Augusta  are  pre¬ 
sented  In  Figures  35,  36,  37*  and  38.  The  data  for  Suf field 
are  the  best  obtained  during  the  project  because  ten  seis¬ 
mographs  were  available.  On  Figure  35,  the  transverse 
ground  motion  Is  plotted  at  Its  correct  azimuth.  The  pres¬ 
ence  of  two  nodes,  with  reversed  polarity  of  either  side.  Is 
quite  clear.  The  Increase  of  high  frequency  energy  in  the 
neighborhood  of  the  nodes  Is  typical.  A  polar  plot  based 
on  the  phase  correlated  by  the  dotted  circle  on  Figure  35 
Is  shown  on  Figure  36.  The  radiation  pattern  for  the  P  wave 
Is  given  on  the  same  figure,  and  the  coincidence  of  a  slight 
bulge  in  the  amplitude  with  the  nodal  lines  for  SH  is  fur¬ 
ther  evidence  that  there  is  a  wave -generating  mechanism 
operative  that  is  superimposed  on  the  basic  compression. 

In  all  cases,  the  radiation  patterns  show  either  two  or  four 
lobes.  In  the  case  of  four-lobed  patterns,  two  of  opposite 
polarity  were  usually  bigger  than  the  other  two. 

The  Florissant  site  is  the  only  one  in  which  enough  shots 
were  fired  to  give  information  on  the  consistency  of  orien¬ 
tations  of  this  radiation  pattern.  A  strong  tendency  was 
found  for  the  nodal  lines  to  be  in  a  fixed  direction  over 
large  areas  of  the  site,  north-northwest  and  south-southeast 
being  the  most  frequently  observed  orientation.  Occasion¬ 
ally  a  pattern  with  a  totally  different  arrangement  of 
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nodes  and  lobes  was  found.  Following  the  initial  discovery 
of  the  SH  radiation  pattern,  the  investigators  expected  to 
find  these  patterns  occurring  at  random.  If  one  postulates 
crack  formation  as  a  primary  source  of  this  motion,  then  a 
tendency  to  a  fixed  orientation  is  plausible.  While  cracking 
from  an  explosion  might  be  expected  to  be  8  somewhat  chaotic 
and  random  process,  it  is  likely  that  in  any  given  test  site 
there  are  directions  of  weakness  in  which  the  material  rup¬ 
tures  most  easily.  In  the  case  of  limestone  or  granite, 
these  directions  might  be  expected  to  be  related  to  the  Joint 
pattern  in  the  region  and  to  regional  tectonic  stresses.  In 
the  case  of  soils,  these  directions  would  be  related  to  the 
structure  of  the  soil  resulting  from  the  process  of  deposi¬ 
tion. 

Radiation  Pattern  as  a  Function  of  Distance. 

Two  further  questions  were  investigated:  Is  the  radiation 
pattern  observed  at  one  distance  maintained  at  greater  dis¬ 
tances,  and  to  what  extent  is  it  possible  to  alter  the  radia¬ 
tion  pattern  by  changing  the  source?  A  definitive  answer 
to  the  first  question  could  have  best  been  obtained  by  sur¬ 
rounding  the  shot  with  at  least  three  rings  of  instruments 
at  different  distances.  Such  a  procedure  would  have  required 
many  more  seismographs  than  were  available,  so  a  compromise 

was  found. 

First,  at  the  Florissant  site  Instruments  were  set  up  at 
two  ranges,  covering  a  quadrant  within  which  a  node  had  been 
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frequently  observed.  The  distances  were  30  and  60  meters 
(S.P.  51  and  52).  Clearly,  If  the  radiation  pattern  Is  not 
maintained  when  the  range  Is  Increased  by  a  factor  of  tv?o, 
the  whole  hypothesis  collapses,  and  further  Investigation 
Is  unwarranted.  The  transverse  component  of  ground  motion 
for  S.P.  51b  Is  presented  In  Figure  39.  It  Is  quite  clear 
that  the  node  In  the  west-northwest  direction  between  Instru¬ 
ments  1  and  2,  3  and  4  Is  expressed  at  both  distances.  Re¬ 
sults  for  52b  are  similar. 

These  data  provide  an  opportunity  to  check  the  method  of 
computing  the  position  of  the  nodal  lines,  viz.,  linear  Inter¬ 
polation  between  positions  of  reversed  polarity,  using  the 
amplitude  of  a  correlated  phase  at  each  station.  For  S.P. 

51b,  if  the  large  peak  on  instrument  6  is  selected  as  the 
phase  to  be  correlated.  Interpolation  between  1  and  2  puts 
the  node  at  24.5°  north  of  west,  while  interpolation  between 
1  and  6  puts  it  20.8°  north  of  west.  If  the  next  turning 
point  is  selected,  that  Is,  the  largest  trough  on  instrument 
6,  the  results  are:  16.5°  using  1  and  2,  and  19.3°  using 
1  and  6. 

If  now  the  records  at  60  meters  are  used,  the  result  cor¬ 
responding  to  the  biggest  peak  on  instrument  5  is  13°  north 
of  west  using  3  and  4,  and  15.4°  using  3  and  5.  The  results 
for  the  following  turning  point  are  27.6°  using  3  and  4, and 
exactly  the  same,  27.6°,  using  3  and  5.  If  all  eight  of 
these  values  are  taken  as  independent  measurements  of  the 
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FIGURE  39.  SH  GROUND  MOTION, 
FLORISSANT,  S.R  51b. 
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angle,  the  mean  value  le  21.2°,  with  a  standard  deviation  of 

4.5°. 

It  may  be  concluded  that  linear  Interpolation  is  Justi¬ 
fied,  since  for  any  one  phase  at  one  distance,  the  result 
for  the  instruments  covering  45°  Is  very  close  to  that  for 
the  instruments  covering  90° •  Most  important  is  the  result 
that  the  mean  position  of  the  node  from  the  data  at  30  meters, 
20.3°  north  of  west  is  very  close  to  the  position  from  the 
data  at  60  meters,  22.2°. 

With  these  satisfactory  results  obtained,  further  exper¬ 
iments  were  conducted  to  determine  if  the  same  pattern  could 
be  observed  as  far  out  as  210  meters.  It  was  not  possible 
to  go  out  this  far  at  the  Florissant  site  without  encounter¬ 
ing  some  pronounced  topographic  features  in  some  directions. 
Because  of  the  shortage  of  available  instruments,  the  data 
from  three  shots  are  combined  in  this  part  of  the  study.  The 
three  shots  were  all  in  the  same  part  of  the  test  site  as 
S.P.  51  and  52,  discussed  above.  Combining  the  results  of 
separate  shots  is  Justified  first  because  the  records  from 
51  and  52  Indicated  not  only  that  the  orientation  of  the 
radiation  pattern  remained  fixed,  but  the  SH  waveforms  were 
fairly  well  reproduced.  Then,  in  order  to  provide  addi¬ 
tional  control,  one  instrument  was  always  kept  in  the  same 
position  relative  to  the  shot,  30  meters  to  the  north.  The 
five  measurements  of  transverse  ground  motion  at  this  posi¬ 
tion  are  presented  in  Figure  40.  Although  there  are  differ- 
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ences  In  detail,  the  major  features,  especially  the  character 
of  the  long  period  Love  wave,  are  reproduced. 

Five  profiles  covering  a  ISO0  arc  from  the  north  through 
west  to  south  could  be  derived  from  the  records  of  S.P.  53b, 
54b,  and  6l.  The  transverse  ground  motion  along  these  lines 
is  given  on  Figures  41  and  42.  Instrument  5  was  to  the  north, 
in  line  with  the  fixed  Instrument  at  30  meters.  Unfortunate¬ 
ly,  It  was  impossible  to  move  out  In  this  direction  more 
than  122  meters.  A  comparison  of  the  records  made  by  Instru¬ 
ment  5  for  S.P.  51b  and  53b,  and  for  5^b  and  61  provides 

additional  evidence  of  reproducibility. 

Timing  control  during  these  experiments  was  not  very  good. 
The  zero- time  device  in  two  of  the  instruments  was  reliable, 
so  that  times  of  selected  events  could  be  taken  from  these 
records.  Then  the  Rayleigh  waves,  v.-hlch,  as  pointed  out  be¬ 
fore,  were  very  much  the  same  In  all  azimuths,  could  be  used 
to  fix  the  time  for  the  other  records  on  which  no  time  break 
appeared.  This  only  worked,  of  course,  for  those  instru¬ 
ments  that  were  at  the  same  distance  from  the  source. 
Therefore,  a  few  of  the  records  for  which  topography  forced 
the  use  of  distances  other  than  the  standard  ones  could  not 

be  timed. 

The  same  radiation  pattern  was  determined  from  these 
three  shots  as  from  the  first  two.  On  some  of  the  records 
the  amplitudes  were  too  small  to  allow  correlations  of 
phases  to  be  made  with  confidence  and  these  were  excluded 
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from  the  analysis.  In  all  cases  the  polarity  on  the  traces 
to  the  south,  southwest,  and  west  was  the  same,  and  reversed 
from  that  to  the  northwest  and  north.  The  mean  position  for 
the  node  was  found  to  be:  at  6l  meters,  24°  north  of  west; 
at  122  meters,  49°  north  of  west,  and  at  213  meters,  34.5° 
north  of  west.  Amplitudes  In  the  southwest  quadrant  did  not 
vary  as  smoothly  as  had  been  noted  previously  at  smaller  dis¬ 
tances,  as  a  very  low  amplitude,  but  with  the  proper  polarity, 
was  observed  on  Instrument  2  at  61  and  122  meters.  Varia¬ 
tions  in  topography  are  most  likely  the  explanation  as  this 
instrument  was  separated  by  a  valley  and  creek  from  the  shot 
point . 

These  data  lead  to  the  final  conclusion  that  the  radia¬ 
tion  patterns  observed  at  short  range  do  propagate  outward. 
This  provides  further  evidence  that  the  polarity  reversals 
are  real,  and  not  apparent  reversals  resulting  from  asymmetri¬ 
cal  dispersion,  and  that  the  Love  waves  are  truly  propagating 
outward  from  the  source  region,  and  are  not  formed  along  the 
path  by  conversion  of  other  wave  types. 

Modification  of  the  Radiation  Patterns. 

Shots  in  Water-Filled  Pits.  Experiments  were  conducted 
at  Florissant  to  test  the  effect  on  the  SH  waves  of  shoot¬ 
ing  in  8 mall  pits  filled  with  water.  These  experiments  must 
be  considered  to  have  failed  because  it  was  impractical  with 
the  facilities  available  to  dig  and  fill  pits  of  sufficient¬ 
ly  great  size  to  provide  a  real  test  of  the  effects.  The 
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original  concept  was  to  shoot  In  a  pit  large  enough  that 
purely  normal  stres.es  within  the  elastic  limit  of  the  sur¬ 
rounding  earth  material,  would  be  applied  to  the  wall,  of 
the  pit.  For  the  small  pits  used,  a  great  deal  of  water 
«as  lifted  vertically,  and  fell  back,  so  that  the  Input  was 
more  like  a  vertically  applied  pulse,  and  the  records  are 
indeed  very  similar  to  those  from  shallow  burled  shots. 

One  series  of  experiment,  will  be  described,  those  at 
S.P.  56.  Shot  56d  was  0.5  pound  charge,  burled  flush  with 
the  surface.  The  crater  produced  by  this  shot  was  dug  out 
to  form  a  hemispherical  pit  with  a  radius  of  0.8  meter.  This 
pit  was  filled  with  water,  and  another  0.5  pound  shot  fired 
at  the  surface  of  the  water.  The  pit  was  then  reshaped  to 
form  an  ellipsoid,  with  the  major  axis  trending  north-south, 
and  twice  as  long  as  the  minor  axis,  refilled,  and  anotner 


surface  shot  fired . 


The  final  test  was  a  shot  in  the  bottom 


of  the  water-filled  ellipsoidal  pit. 

The  analysis  of  the  data  is  summarized  In  Figures  43  to 
46.  The  efrect  on  all  components  of  shooting  In  the  water 
was  minor.  The  shot  at  the  bottom  of  the  pit  produced 
records  very  similar  to  the  original  tamped  shot,  but  with 
greatly  Increased  amplitudes  on  all  component. .  In  particu- 
lar  the  3H  motion  from  this  shot  was  very  much  like  the 
original  shot,  so  tl»t  repeated  loading  of  the  same  source 
region  results  In  reproducible  motion.  This  fact  must  be 
reconciled  with  the  hypothesis  that  crack  formation  is 
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FIGURE  43 

RADIATION  OF  SURFACE  WAVES 


S.P  56  d 
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FIGURE  45 
ALTERATION  OF  SH 
BY  SHOOTING  IN  THE  PITS 
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responsible  for  the  SH  motion.  Apparently  once  an  under¬ 
ground  crack  pattern  has  been  established,  repeated  load¬ 
ing  results  in  further  motion  along  the  same  cracks.  The 
fluid  in  the  pit  fills  the  crack3,  and  the  pressure  from 
the  shot  is  transmitted  through  the  fluid  to  produce  fur¬ 
ther  motion  of  the  crack  walls,  as  well  as  further  extension 
of  the  crack. 

Shooting  at  the  top  of  the  water-filled  pit  reduces 
the  amplitude  of  all  events  in  the  same  manner  as  a  poorly 
tamped  surface  shot  in  soil.  The  SH  motion  is  altered  in 
some  directions,  and  the  effect  on  SH  is  greater  than  on 
the  Rayleigh  wave,  though,  it  must  be  emphasized,  the  ob¬ 
served  effects  are  minor.  Even  the  minor  effects  are  diffi¬ 
cult  to  interpret,  since  possibly  effects  of  hole  fatigue 
are  included.  The  Rayleigh  wave  is  used  as  a  basis  of  com¬ 
parison  because  its  waveform  is  so  well  known  at  this  site, 
it  is  symmetric  in  its  distribution,  and  it  is  reproducible 
from  shot  to  shot.  Shot  56d  was  used  as  the  calibration 
shot.  A  phase  in  the  vertical  component  of  the  Rayleigh 
wave  (Rz),  the  longitudinal  component  of  the  Rayleigh  wave 
(Rl),  and  in  the  Love  wave  (SH)  that  could  be  correlated  on 
all  of  the  seismograms  was  selected  and  its  amplitude  read. 
The  radiation  of  these  three  phases  for  the  calibration  shot 
is  shown  in  Figure  43.  The  two-lobed  SH  pattern,  with  the 
typical  northwest-southnast  orientation  of  the  nodes,  and 
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the  almost  circular  pattern  for  the  two  Rayleigh  components 
is  seen. 

The  pattern  for  56e,  the  top  center  of  the  circular 
pit,  now  shown  here,  was  very  similar.  The  minimum  in 
to  the  southwest  is  more  pronounced.  The  shot  in  the  top 
of  the  elliptical  pit,  56f,  produced  a  poorly  defined  four- 
lobed  SH  pattern.  In  particular,  the  reading  on  the  north¬ 
east  instrument  on  which  the  Interpretation  as  a  four-lobed, 
rather  than  two-lobed,  pattern  depends  is  very  doubtful 
because  of  a  small  amplitude  and  great  difficulty  in  corre¬ 
lating  the  selected  phase.  When  the  shot  was  placed  at  the 
bottom  of  this  same  pit,  a  distorted  two-lobed  pattern  was 
obtained.  The  SH  amplitude  on  the  north  instrument  (Num¬ 
ber  l)  in  shot  56d  was  taken  as  the  unit  amplitude,  and  the 
patterns  for  all  four  shots,  referred  to  this  unit,  are 
shown  in  Figure  44. 

The  alteration  of  the  correlated  phase  relative  to  the 
original  shot  was  measured  by  dividing  the  amplitude  on  each 
instrument  by  the  value  on  the  same  instrument  for  56d.  If 
there  were  no  effect  at  all,  the  result  would  be  a  circle 
of  unit  radius,  and  if  the  radiation  were  unchanged  but  the 
amplitudes  were  altered,  the  plot  would  be  a  circle  with 
appropriate  radius.  As  seen  in  Figure  45,  the  actual  result 
for  56e  is  a  slight  increase  in  amplitudes  in  the  northwest 
direction  and  a  decrease  to  the  northeast  and  southwest. 

The  bottom  shot,  56h,  produced  SH  motion  roughly  twice  as 
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larse  as  the  original  shot,  but  the  signal  to  the  southeast 
was  three  times  the  original,  while  that  to  the  northwest 
was  only  about  one  and  one-half  times  the  original,  and  had 
the  opposite  sign.  The  reversal  In  sign  at  a  station  near 
a  node  Is  not  surprising,  and  Is  also  observed  on  56f. 

The  three  shots  taken  together  Indicate  a  tendency  for 
a  relative  strengthening  of  the  signal  In  the  northwest- 
southeast  direction,  but  this  tendency  has  no  obvious  con¬ 
nection  with  the  pit,  which  was  either  circular,  or  elong¬ 
ated  north  and  south.  It  Is  recommended  that  further  work 
In  which  the  shots  are  fired  In  bodies  of  water  big  enough 
that  any  directly  produced  shear  motion  Is  decoupled  from 
the  earth  be  carried  out. 

The  amplitude  of  the  Rayleigh  waves  were  altered  In 
roughly  the  same  manner  as  the  Love  waves,  that  Is,  those 
shots  that  produced  smaller  Love  waves  also  produced  smaller 
Rayleigh  waves,  though,  of  course,  the  Rayleigh  wave  radia¬ 
tion  remained  quite  symmetric.  The  Rayleigh  wave  radiation 
normalized  to  the  north  Instrument  56d,  Is  shown  In  Figure 
46,  in  which  It  Is  seen  that  the  bottom  shot  produced  Ray¬ 
leigh  waves  with  twice  the  amplitude  of  the  original  burled 
shot,  while  the  shots  at  the  top  of  the  pits  produced  appre¬ 
ciably  smaller  motion.  No  evidence  of  any  effect  of  the 
shape  of  the  pit  on  the  radiation  pattern  can  be  seen.  ■ 

A  final  comparison  that  can  be  made  Is  the  enhance¬ 
ment  of  SH  motion  relative  to  the  enhancement  of  the  Ray- 
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leigh  wave.  To  measure  this  the  ratio  of  SH  for  the  pit 
shots  divided  by  the  values  at  the  same  Instruments  for  the 
calibration  shot  (the  values  on  Figure  44)  were  divided  by 
the  corresponding  ratios  for  the  Rayleigh  wave  (not  plotted 
here).  The  values  are  scattered  somewhat,  but  on  the  aver¬ 
age,  the  value  for  the  top  of  the  circular  pit  was  1.91*  for 
the  top  of  the  elliptical  pit  2.03*  and  for  the  bottom  of 
the  elliptical  pit,  0.95.  Therefore,  rather  than  reducing 
the  production  of  SH  motion,  as  had  been  anticipated,  the 
shots  in  the  water  actually  doubled  it  relative  to  the  Ray¬ 
leigh  wave,  while  the  shot  in  the  bottom  of  the  water  yielded 
the  same  proportions  of  Rayleigh  and  SH  motion. 

SH  Motion  from  Ripple -Fired  Shots.  The  Rayleigh  wave 
motion  produced  by  the  ripple-fired,  multi-hole  shots  was 
discussed  in  Chapter  IV.  The  SH  motion  generated  by  the 
horizontal  patterns  will  be  described  here.  The  experiments 
were  outlined  in  the  previous  chapter.  The  method  that  has 
been  used  to  study  SH  radiation  does  not  work  as  well  for 
these  seismograms  because  the  superposition  of  the  signals 
from  the  component  parts  of  the  charge  produces  waveforms 
that  cannot  be  readily  correlated  around  the  circle  of  in¬ 
struments.  As  explained  in  the  discussion  of  Rayleigh  waves, 
since  the  signal  from  a  single  shot  is  almost  identical  in 
all  azimuths,  the  superposition  of  this  waveform  with  proper 
time  shifts  will  explain  the  recorded  motion.  Because,  as 
has  been  seen  above,  the  SH  radiation  from  a  single  shot  is 
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fundamentally  complicated  and  largely  unpredictable,  analyti¬ 
cal  euperpoeltlon  la  Impractical,  even  though  this  le  un¬ 
doubtedly  the  way  to  explain  the  obeerved  motion. 

The  calibration  shots  for  the  rlpple-flred  series,  S.P. 

69a  and  b,  produced  a  four-lobed  pattern  (Figure  47)  almost 
Identical  to  that  predicted  by  a  formula  given  by  Knopoff 
and  Gilbert  for  the  first  motion  due  to  a  crack  growing  rrcm 
S70°W  to  N  70°  E  at  twice  the  shear  velocity  (see  Figure  12 
in  Paper  I).  The  two-hole  shot  fired  with  the  shorter  time 
delay  (equivalent  to  propagation  velocity  of  190  meters/sec. , 
S.P. 70)  produced  smaller  SH  motion  on  all  components.  The 
character  of  the  motion  In  the  direction  away  from  which  the 
shot  was  propagating  was  completely  altered,  and  no  corre¬ 
lation  with  previous  records  or  the  other  records  for  this 
shot  could  be  made.  The  phase  that  had  been  correlated  In 
analyzing  S.P.  69  could  be  picked  out  on  the  records  toward 
which  the  shot  propagated,  with  the  same  polarity.  On  the 
shot  fired  with  slow  velocity  of  91  meters/sec.  (S.P.71), 
the  character  was  much  the  same  as  for  the  shorter  delays, 
but  It  was  possible  to  carry  a  correlation  through.  The 
resulting  pattern  of  radiation  (Figure  47)  consisted  of  only 
two  lobes,  a  small  one  which  covered  the  two  smaller  lobes 
from  69a,  and  a  large  one  covering  roughly  the  two  larger 
lobes  from  that  calibration  shot.  Nodal  lines  to  the  north 

and  to  the  southeast  were  preserved. 

Tracings  of  the  original  seismograms  from  three  five-hole 
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shots,  and  from  the  tw  eight-hole,  two-llne  pattern  are 
shown  in  Figure  48.  For  all  shots  the  Instruments  were 
placed  in  a  circle  with  Number  1  to  the  north  (roughly), and 
60°  spacing,  with  2  to  the  northwest,  3  southwest,  4  south, 
etc.  For  the  five-hole  patterns,  the  charges,  spaced  at 
five  foot  intervals,  were  aligned  north  to  south,  and  det¬ 
onated  from  south  to  north  (from  Instrument  4  to  Instrument 
1) .  Geometrically,  2  and  6  were  symmetrically  placed  with 
respect  to  the  line,  as  were  3  and  5.  For  the  "simple 
couple"  type  shots  (S.P.  72a  and  b),  the  eight  shots  placed 
on  two  parallel  lines,  five  feet  apart,  extending  north  to 
south.  Detonation  was  Initiated  simultaneously  at  the  south 
end  of  the  west  line  and  the  north  end  of  the  east  line.  In 
examining  the  seismograms  In  Figure  48,  It  must  be  recalled 
that  the  static  magnification  of  Instruments  3,  4  and  5  la 
about  1.8  times  that  of  the  other  three. 

The  shot  in  which  all  the  charges  were  fired  simultan¬ 
eously  (67a)  gives  the  basic  pattern  for  the  five-hole  shots. 
It  may  be  seen  that  the  motion  recorded  in  line  with  the 
pattern  (at  1  and  4)  is  a  minimum.  The  Love  wave  pulse 
can  be  correlated  on  all  records  exoept  1.  The  records  at 
the  symmetric  points  2  and  6  are  very  much  the  same,  though 
6  has  larger  amplitudes,  and  the  other  symmetric  pair,  3 
and  5,  have  easily  correlated  phases  of  similar  amplitudes 
but  with  reversed  polarity.  When  the  same  line  of  shots  is 
fired  at  an  effective  propagation  velocity  of  290  meters/sec. 
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FIGURE  48 


SH  MOTION  FROM  HORIZONTAL  PATTERNS 


-127- 


(67b)  the  same  waveforms  are  observed,  only  a  bit  more  ex¬ 
tended  in  time  as  a  result  of  the  duration  of  the  input.  No 
significant  change  in  radiation  is  noted.  The  slower  speed, 
110  meters/sec.  (68a),  resulted  in  marked  reduction  of  the 
large  amplitudes  at  instruments  3  and  5,  and  proportionally 
greater  motion  on  the  other  instruments;  in  other  words, 
more  uniform  radiation. 

Because  of  instrument  difficulties,  the  simple  couple 
experiments  were  incompletely  recorded.  The  most  remark¬ 
able  effect  observed  is  the  increase  in  amplitude  in-line 
with  the  pattern.  In  particular,  the  seismogram  at  the 
south  end.  Number  4,  is  very  large.  When  it  is  recalled 
that  this  was  the  direction  of  minimum  motion  for  the  single 
line  pattern,  the  enhancement  must  be  regarded  as  an  effect 
of  the  pattern.  However,  if  the  pattern  truly  simulated  a 
couple,  this  should  be  the  direction  of  a  node  of  SH  motion. 
This  result  is  unexplained.  The  extended  duration  of  the 
motion  is  also  to  be  noted. 


CHAPTER  VI 


THE  EFFECT  OF  SOURCE  DEPTH  AND  SHOT  POINT 
MEDIUM  ON  THE  SEISMIC  SIGNAL 


Knowledge  of  the  properties  of  explosion-generated 
seismic  waves  Is  largely  empirical  based  on  the  analysis  of 
countless  seismograms  produced  by  explosions  covering  a 
tremendous  range  of  yields  and  fired  under  a  wide  variety 
of  conditions.  Attempts  to  understand  this  large  body  of 
data  have  led  to  theoretical  developments  concerning  the 
processes  of  wave  formation  and  wave  propagation  that  are 
the  foundation  of  this  branch  of  seismology.  Many  aspects 
of  this  theory  are  well-developed,  byt  because  of  the  com¬ 
plex  nature  of  the  materials  of  which  the  earth  Is  conposed 
and  the  complex  structure  of  the  medium  through  which  the 
waves  must  propagate  to  the  seismograph,  the  Interpretation 
of  particular  observations  Is  often  difficult.  In  particu¬ 
lar,  the  seismologist  Is  called  upon  to  Isolate  the  effects 
of  a  number  of  variable  parameters  that  are  part  of  the  cir¬ 
cumstances  surrounding  any  test  explosion. 

Data  concerning  the  effects  of  several  of  these  vari¬ 
ables,  depth  of  source,  medium  In  which  the  source  is  lo¬ 
cated,  and  charge  size,  were  gathered  during  this  project. 
The  results  on  the  effect  of  charge  size  are  of  doubtful 
value  for  two  reasons:  the  range  of  yields  available  was 
not  sufficiently  great  to  adequately  define  significant 
trends;  and  during  some  of  the  tests  that  were  specifically 
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deslgned  to  give  data  on  this  question  the  experimenters 
were  unfortunate  in  getting  some  bad  dynamite,  which  did  not 
conqpletely  detonate.  Some  references  to  the  effect  of  charge 
size  are  Included  In  the  other  sections  of  this  report,  but 
no  systematic  evaluation  is  possible.  The  conclusions  con¬ 
cerning  the  other  two  parameters  will  be  discussed  In  thia 
chapter.  For  a  summary  of  the  present  state  of  knowledge  of 
the  factors  affecting  the  close-in  properties  of  the  seismic 
signal,  see  Klssllnger  (1963). 

Effect  of  Source  Depth  on  the  Seismic  Signal. 

The  depth  of  the  charge  below  the  surface  will  affect 
both  body  and  surface  waves.  Surface  waves  depend  for  their 
existence  on  the  presence  of  a  free  surface,  and,  in  general, 
their  amplitude  decreases  as  the  source  depth  Increases. 

The  effect  of  Increased  source  depth  is  not  as  simple  for 
body  waves.  If  the  shot  point  medium  is  homogeneous,  the 
effect  of  Increased  shot  depth  is  to  reduce  the  amplitude 
of  the  direct  P-wave,  This  results  from  increased  ambient 
stresses  at  the  shot  which  make  the  cavity  size  smaller. 

This  decrease  in  cavity  size  also  causes  an  increase  in  the 
frequency  of  the  direct  P-wave.  In  a  real  experimental 
situation,  the  properties  of  the  medium  at  the  shot  will 
change  with  depth,  and  thlB  will  override  the  effect  of  in¬ 
creased  stress.  That  is,  for  a  limited  range  of  depths, 

P-wave  amplitude  and  frequency  will  depend  more  on  the 
medium  than  on  the  depth. 
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Effect  on  P-Waves,  The  predominant  frequency  In  the 
p-wave  has  been  found  to  be  affected  markedly  by  the  medium 
at  the  depth  of  the  charge.  The  analysis  of  a  step  function 
In  pressure  applied  to  the  walls  of  a  spherical  cavity  in  an 
Infinite  medium  Indicates  that  at  distances  beyond  two  or 
three  cavity  radii,  the  P-wave  should  be  damped  sinusoid 
with  angular  frequency  given  by 


6J  «=(q  +  q^)^  V 

(1  +  2qr  a 

wherei  V  Is  the  P-wave  velocity;  a  Is  the  cavity  radius; 
q  o  t  the  ratio  of  shear  modulus  to  Lame's  constant 

(Sharpe,  1942;  Blake,  1952,  Kisslinger,  1963a).  The  damp¬ 
ing  is  given  by 


1  +  2, 


where  S  is  the  fraction  of  critical  damping.  For  large 
values  of  Poisson's  ratio  (small  q)  low  frequencies  and 
small  damping  are  to  be  expected.  The  radius,  a,  of  the 
equivalent  cavity  Is  usually  not  known  for  a  given  medium 
and  yield. 

At  Florissant  a  series  of  shots  In  two  holes,  one  14,6 
meters  deep  and  the  other  38. 1  meters  deep,  were  fired  to 
gather  data  for  P-  and  S-wave  velocities  and  the  study  of 
the  generation  of  Rayleigh  waves.  These  shots  also  yielded 
data  useful  for  checking  the  theory  of  the  waveform  of  the 


direct  P-wave. 
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At  this  test  site  the  velocity  of  the  P-wave  in  the 
clay-loess  overlying  the  Mississipplan  limestone  changes  at 
the  water  table  whereas  the  S-wave  velocity  does  not  show  a 
marked  change  at  this  boundary.  This  results  in  the  high 
values  of  Poisson's  ratio  presented  in  Chapter  III.  Above 
the  water  table  the  P-wave  velocity  has  an  average  value  of 
460  meters  per  second.  For  this  part  of  the  section  A  = 

3.3 /a.  ,  consequently  the  frequency  of  the  P-wave  is  f  **co/2W 
=  57.5/a  cps,  where  a  is  in  meters,  jf  is  0.43,  and  the  damp¬ 
ing  ratio  (peak  to  trough)  is  4.53.  Below  the  water  table 
the  P-wave  velocity  is  1800  meters  per  second,  /  =  42.5^m  , 

f  =  84.9/a  cps,  5  is  0.15,  and  the  damping  ratio  is  1.61, 

These  predicted  values  may  be  compared  with  those  observed, 
presented  in  Table  4.  These  data  were  taken  from  records  at 
a  range  of  2.5  meters  from  1.5  pound  shots  in  a  single  hole 
at  depths  from  38  to  7.5  meterB .  For  the  shots  deeper  than 
l4  meters  both  frequencies  and  damping  ratios  are  close  to 
those  predicted  for  a  cavity  radius  of  one  meter.  It  is 
assumed  that  the  damping  associated  with  the  medium  at  the 
source  is  preserved  in  the  signal  observed  at  the  surface. 

The  same  radius  gives  good  agreement  for  frequencies,  but  not 

for  the  damping  ratios  for  the  shallower  shots.  The  vel¬ 

ocities  used  for  the  calculations  are  averages  for  the  layer, 
and  there  may  be  unknown  variations  in  the  velocities  in 
thin  layers  at  any  level. 
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Table  4 


FREQUENCY  AND  DAMPINO  OF  P-WAVE 
FLORISSANT  SITE,  RECORDED  AT  2.5  METERS 


Depth 

(meters) 


Damping  Ratio 
Frequency  (First  Peak 

( cps )  to  Trough) 


38  77-80  1.30 

32  91  1.61 

26  77  1.43 

20  59 

14  77  0.96 

8  59  1.50 


The  observed  P-wave  motion  over  the  shot  point  from  a 
series  of  0.5  pound  charges  at  various  depths  In  the  14.6 
meter  hole  showed  considerable  variation  in  waveform.  The 
damping  ratio  was  much  less  than  that  predicted  for  moBt  of 
the  shots.  This  was  thought  to  be  due  to  the  interference 
of  the  direct  wave  and  the  primary  reflections  from  the 
known  boundaries  in  the  clay  section.  To  check  this,  theo¬ 
retical  seismograms  of  the  P-motion  were  made.  The  waveforms 
used  as  input  for  these  theoretical  seismograms  were  constant- 
frequency  damped  sinusoids  with  the  first  peak  equal  to 
unity.  The  theoretical  frequency  and  damping  given  above 
were  used.  The  geologic  model  used  for  the  calculations, 
given  in  the  table  below,  was  made  from  observed  uphole 
times  from  the  deep  holes. 

Thickness  P-Velocity 
(meters)  (meters/sec . ) 


6.1 

9.2 

30.5 

OO 


305 

518 

1830 

3100 
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In  calculating  the  theoretical  seismograms,  transmission 
and  reflection  coefficients  at  each  interface  were  calcu¬ 
lated,  assuming  normal  incidence.  All  primary  reflections 
from  boundaries  below  the  shot  were  added  to  the  directly 
transmitted  wave  at  the  proper  time  delay. 

The  comparison  of  the  observed  and  calculated  motions 
from  tnree  of  the  shots  in  the  third  layer  are  shown  in 
Figure  49.  The  general  agreement  is  good.  The  theoretical 
motion  for  the  source  at  26  meters  shows  stronger  motion 
after  the  arrival  time  of  the  reflection  from  the  top  of  the 
Missslssippian  limestone  than  is  seen  on  the  field  record. 

Theoretical  P  motion  for  shots  in  the  upper  layers  were 
calculated  using  a  waveform  with  a  damping  ratio  of  4.53.  The 
theoretical  motion  did  not  resemble  the  observed  motion  be¬ 
cause  the  damping  was  much  too  large.  Another  set  was  cal¬ 
culated  using  a  waveform  with  the  theoretical  frequency  of 

57.5  cps,  but  the  damping  ratio  equal  to  1.6l.  Figure  ^0 
shows  the  comparison  of  the  observed  and  calculated  for 
shots  in  the  middle  layer.  The  agreement  for  the  shots  at 

14.6  and  11.6  meters  is  good.  There  is  a  zero  shift  on  the 
observed  traces  resulting  from  the  strong  motion  hitting  the 
seismometer.  There  is  fair  agreement  for  the  shot  at  9.3- 
meters  but  the  agreement  is  very  poor  for  the  shot  at  7.6 
meters.  Peaks  and  troughs  are  observed  at  approximately 
the  correct  times,  but  the  relative  amplitudes  are  much 
reduced.  It  is  evident  that  the  expected  reflections  are 
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not  well  developed  for  the  shallow  shot.  Whether  this  ab¬ 
sence  of  strong  down-going  energy  from  the  shot  is  due  to 
losses  In  the  shallow  zone  or  scattering  from  the  cavities 
left  by  the  prior  deeper  shots  Is  not  known.  The  agreement 
Indicates  that  the  damping  used  In  the  calculations  Is  more 
nearly  correct  than  the  value  predicted  from  theory.  The 
combined  effect  of  decreased  damping  and  consideration  of 
the  primary  reflected  eventB  along  with  the  direct  wave  give 
a  reasonable  explanation  for  the  oscillatory  nature  of  the 
P-motion  at  these  Intermediate  depths. 

The  comparison  of  the  theoretical  Belsmograms  and  the 
observed  motions  for  shots  In  the  upper  layer  are  shown  In 
Figure  51.  The  agreement  Is  not  as  good  as  at  the  greater 
depths.  The  most  striking  difference  1b  the  apparent  ab¬ 
sence  of  a  reflection  from  the  boundary  at  a  depth  of  15.3 
meters  (water  table)  on  the  observed  trace.  ThlB  reflection 
is  strong  on  the  calculated  traces.  The  observed  frequen¬ 
cies  are  lower  than  those  in  the  calculated  motion. 

The  general  agreement  of  the  theoretical  seismograms  and 
the  observed  traces  from  shots  in  the  two  lower  layers  can 
be  interpreted  to  mean  that  the  P-wave  from  these  contained 
explosions  haB  the  waveform  of  a  damped  sinusoid  with  a 
frequency  predicted  by  the  theory  discussed  at  the  beginning 
of  this  section.  However,  sb  the  source  moved  closer  to 
the  surface,  lower  values  of  damping  than  that  given  by  the 
theory  were  required  to  get  reasonable  agreement.  The  pro- 
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longed  oscillatory  nature  of  the  P-motlon  observed  from 
shots  at  certain  depths  Is  adequately  explained  by  the  com¬ 
bined  effect  of  the  primary  reflections  from  the  near  sur¬ 
face  boundaries  and  the  oscillatory  nature  of  the  initial 
waveform. 

At  the  very  end  of  the  research  program,  a  series  of 
larger  charges  (4.5  pounds)  was  fired  under  conditions  sim¬ 
ilar  to  the  deep  shots  discussed  above .  The  vertical  com¬ 
ponent  of  ground  motion  over  the  shot  point  from  4.5  pounds 
at  38»  32.2  and  19.9  meters  is  shown  in  Figure  52.  The 
arrival  time  of  the  reflection  from  the  top  of  the  lime¬ 
stone  is  also  shown.  The  waveforms  are  not  oscillatory  as 
in  the  previous  examples,  but  appear  as  highly  damped  pulses. 
The  damping  of  the  waveform  should  depend  only  upon  the 
elastic  properties  of  the  medium  and  not  upon  the  size  of 
the  charge.  The  observed  decrease  in  frequency  is  expected 
with  the  increase  in  cavity  size,  but  the  decrease  to  8  cps 
from  85  cps  for  a  charge  increase  of  three  times  is  too 
great  to  be  explained  on  any  such  simple  basis.  These  data 
are  completely  anomalous  for  the  test  site.  Admittedly  it 
proved  discouraging,  after  the  excellent  results  for  the 
first  series  discussed,  to  encounter  a  set  of  data  thf.t  did 
not  fit  the  same  simple  theory.  It  is  evident  that  a  fur¬ 
ther  study  of  the  possible  effects  of  yield  on  damping  as 
well  as  frequency  of  the  primary  signal  is  required. 

The  P-wave  amplitudes  were  studied  from  one  sequence 


ARRIVAL  TIME  OF  REFLECTION 
FROM  TOP  MISSI  SSI  PPI  AN 


I 


OBSERVED  P  MOTION  OVER  THE  SOURCE 

4.5  lb.  CHARGES 
-FLORISSANT- 


-140- 


of  charges.  The  amplitudes  measured  (vertical  component) 
directly  above  a  sequence  of  0.5  pound  shots  fired  at  depth 
of  14.5  to  1.2  meters  closely  fit  a  decay  law  of  d*1*5. 

The  effect  of  distance  is  combined  with  depth,  and  undoubted¬ 
ly  is  predominant.  This  is  somewhat  more  rapid  attenuation 
than  ordinary  spherical  wave  propagation  implies,  but  a  more 
rapid  rate  than  r”1  for  the  peak  displacement  is  expected 
even  for  perfectly  elastic  behavior  in  the  region  near  the 
source,  as  shown  by  J.  Vanek  (1959).  Thus  the  effects  of 
viscous  losses,  transitional  propagation,  and  depth  are  com¬ 
bined,  masking  the  effect  of  any  one  of  these. 

Effect  on  Rayleigh  Waves.  The  effect  of  depth  on  Ray¬ 
leigh  waves  is  well  illustrated  in  a  sequence  of  three  shots 
at  depths  of  0.82,  2.72  and  5.3  meters  recorded  at  25  meters 
at  Florissant  (S.P.  12).  This  motion,  shown  in  Figure  53, 
shows  the  decrease  in  amplitude  with  increasing  shot  depth 
even  though  the  charge  size  of  the  deepest  shot  is  twice 
that  of  the  shallow  shot.  Another  interesting  feature  of 
these  shots  is  the  emergence,  as  the  depth  Increases,  of  an 
event  preceding  the  Rayleigh  motion.  This  event  is  most 
probably  a  higher  mode,  and  the  deepest  shot  may  have  been 
located  near  an  antinode  for  this  mode .  The  presence  of 
the  Mgi  mode  at  this  site  is  discussed  in  Chapter  III. 

In  an  effort  to  gain  a  clearer  understanding  of  the 
generation  of  Rayleigh  waves  from  a  contained  explosion, 
several  sequences  of  shots  in  deep  holes  were  fired  at  the 
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Florissant  test  site.  One  sequence  of  0.5  pound  charges, 
at  various  depths  between  1.0  and  14.6  meters,  was  Instru¬ 
mented  with  three-component  Sprengnether  seismographs 
aligned  along  a  profile  at  0,  2.5,  5,  7.5,  10,  15  and  25 
meters.  This  was  supplemented  with  the  SIE  refraction 
equipment.  Two  vertical  geophones  were  placed  beside  each 
Sprengnether  between  0  and  15  meters.  The  output  of  one 
geophone  at  each  distance  was  passed  through  a  0  to  13  cps 
filter  and  the  other  through  14-100  cps.  This  was  done  In 
an  effort  to  record  on  the  high  gain  equipment  any  long 
period  motion  that  might  be  present  In  the  source  region. 

The  vertical  and  radial  components  of  ground  motion  at 
2.5  meters  are  shown  in  Figure  54,  and  the  motion  at  15 
meters  In  Figure  55.  At  2.5  meters  the  motion  from  the  3 
meter  shot  Is  dominated  by  the  Rayleigh  wave  but  this  event 
rapidly  disappears  as  the  shot  depth  Increases.  At  15  meters 
the  motion  from  the  shallow  shot  is  largely  Rayleigh  wave. 

The  vertical  component  shows  the  rapid  decay  of  amplitude  of 
this  event  which  appears  to  be  totally  absent  at  this  dis¬ 
tance  from  the  shots  deeper  than  7.6  meters. 

The  motion  from  these  shots  was  Fourier  analyzed  to 
obtain  the  Bpectrum  at  one  distance  for  all  depths  and  at 
all  diBtanceB  for  one  depth.  The  result  of  these  analyses 
are  shown  In  Figures  56  and  57.  Figure  56  contains  the 
spectra  of  the  vertical  and  radial  components  for  various 
distances  from  0.5  pounds  at  3  meters.  The  peak  between 
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10  and  20  cps  corresponds  to  the  Rayleigh  wave.  The  peak  at 
37  cps  on  the  vertical  component  at  2,5  meters  from  the  shot 
point  Is  from  the  P-wave.  As  distance  Increases  the  Ray¬ 
leigh  energy  dominates  the  spectrum  of  the  vertical  compo¬ 
nents.  The  sharp  peak  at  approximately  40  cps  on  the  radial 
component  may  be  associated  with  a  higher  mode  of  the  Ray¬ 
leigh  wave.  This  peak  is  prominent  on  the  radial  component. 
Figure  57,  at  7.5  meters  for  the  shot  at  4.6  meters  but  is 
absent  from  the  shot  at  6.1  meters  and  begins  to  reappear  on 
the  deeper  shots.  This  would  mean  that  the  shots  at  3  and 
4.6  meters  are  near  an  antinode  and  the  shot  at  6.1  meters 
is  near  a  node.  Thir  figure  displays  strikingly  the  change 
in  the  energy  in  the  Rayleigh  wave  as  depth  increases.  The 
Rayleigh  peak  diminishes  to  the  point  that  for  shots  at  7*6 
meters  and  deeper  the  spectral  amplitude  is  in  the  noise 

level  in  the  range  of  10  to  20  cps. 

These  data  support  the  result  derived  by  Nakano  (1925), 
that  there  is  a  zone  on  the  free  surface,  the  radius  of 
which  depends  upon  the  depth  of  the  source  and  the  veloci¬ 
ties  of  the  medium,  within  which  the  Rayleigh  wave  does  not 
exist.  This  distance  is  given  by  (Ewing,  Jard^zky  and 
Press,  1957) 

x 

where:  h  is  the  depth  of  the  shot, 

cr  is  the  Rayleigh  wave  velocity,  and 

e<  is  the  P-wave  velocity. 
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Under  the  assumption  that  the  Rayleigh  wave  velocity  Is 
approximately  equal  to  the  shear  velocity  of  the  upper  lay¬ 
er,  a  good  approximation  for  12  cps  (see  Florissant  phase 
velocity  curve.  Figure  7)»  the  critical  distances  for  the 
shots  In  this  sequence  are: 

Depth  Distance 

(meters)  (meters) 

3.05  1.^7 

4.60  2.22 

6.10  2.95 

7.60  3.67 

9.10  4.4o 

11.60  5.60 

14.60  7.05 

It  has  been  shown  that  the  Rayleigh  energy  Is  very  much  re¬ 
duced  on  the  records  at  7.5  meters  from  the  shots  deeper 
than  6.1  meters,  yet  the  table  showB  that  this  distance  Is 
near  the  critical  distance  for  the  14.6  meter  shot  but  well 
beyond  the  critical  value  for  all  others.  Pekerls  and 
Llfson  (1957)  calculated  curves  of  vertical  and  horizontal 
displacement  at  the  free  surface  sb  a  function  of  the  ratio 
of  radial  distance  to  depth  of  source  for  a  point  source. 
These  curves  showed  that  the  Rayleigh  wave  was  not  identi¬ 
fiable  for  distances  less  than  or  equal  to  depth  of  source 
divided  by  the  square  root  of  two  and  was  barely  visible 
on  the  vertical  component  at  five  times  the  depth.  Thus 
for  a  shot  7.6  meters  deep,  the  vertical  component  of  the 
Rayleigh  wave  would  be  poorly  developed  at  38  meters  from 
the  source  and  probably  would  not  be  recorded  on  the  low 
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magnification  Sprengnethers  used  in  this  experiment. 

However,  the  filtered  high  gain  records  obtained  with 
the  SIE  refraction  equipment  contained  long-period  motion, 
approximately  12  cps,  even  over  the  source  for  all  depths  of 
shots  betwen  1  and  14.6  meters.  This  is  shown  for  the  3 
meter  shot  in  Figure  58.  This  long  period  motion  is  seen 
to  step  out  across  the  spread.  The  possibility  that  this 
long  period  motion  was  actually  the  response  of  the  filters 
and  amplifiers  to  an  impulsive  input  was  considered  and  re¬ 
jected  for  the  following  reasons: 

l)  the  same  motion  is  seen  on  the  low  sensitivity  traces 
produced  by  the  completely  mechanical  Sprengnether  portable 
seismographs;  2)  the  impulse  response  of  the  recording  sys¬ 
tem  is  known  (shown  in  the  figure),  and  while  there  is  a 
fairly  long  transient,  for  filter  setting  0-13,  it  1b  not 
oscillatory  and  does  not  resemble  the  recorded  motion. 

These  records  indicate  that  even  inside  the  zone  in  which 
the  Rayleigh  wave  should  not  exist,  there  is  oscillatory 
motion  of  long  duration.  Presumably  thlB  will  be  true  as 
long  as  the  charge  is  buried  deeply  enough  that  no  cratering 
takes  place. 

If  it  la  tr"e  that  oscillatory  motion  of  the  free  surface 
does  exist  at  very  short  range,  the  characteristics  of  this 
motion  need  to  be  explained.  Whereas  the  oscillatory  nature 
of  the  P-wave  1b  predicted  by  simple  theory,  and  the  fre¬ 
quency  of  this  wave  is  known  in  terms  of  the  properties  of 
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the  medium,  no  similar  theory  for  the  Rayleigh  wave  at  very 
short  range  is  known. 

Dix  (1955)  discussed  a  possible  mechanism  for  the  gen¬ 
eration  of  "long  waves"  from  explosions.  His  hypothesis 
is  as  follows:  a)  the  spreading  spherical  pulses  from  a 
buried  compressional  source  are  converted  into  non-spherical 
secondary  disturbances  on  striking  the  free  surface  (essen¬ 
tially  a  diffraction  phenomenon);  b)  because  of  the  non- 
spherical  character,  the  propagation  is  no  longer  simply  in 
the  outward  direction  but  is  generally  in  every  available 
direction,  therefore  an  appreciable  amount  of  energy  goes 
back  toward  the  source  or  toward  a  vertical  axis  through  the 
source;  c)  thus,  energy  is  held  near  the  source  and  near 
the  free  surface  for  a  long  time  compared  with  the  time  it 
takes  for  energy  of  the  signal  to  be  transmitted  away  from 
the  source  neighborhood;  and  d)  this  long  held  energy 
spreads  outward  in  every  available  mode,  one  of  which  is 
the  surface  Rayleigh  wave.  Kobayashi  and  Takeuchi  (1957) 
considered  surface  motion  from  a  buried  compressional  line 
source.  They  calculated  profiles  of  surface  displacements 
at  certain  fixed  times  from  an  impulsive  source.  They 
found  that  the  characteristic  surface  profile,  one  which 
persists  as  the  basic  waveform  for  all  time,  is  formed  be¬ 
tween  a  time  equal  to  twice  the  depth  divided  by  the  P- 

wave  velocity  and  five  times  the  same  ratio.  They  related 
the  surface  motion  before  and  after  this  time  to  the 
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mechanism  proposed  by  Dix.  The  presence  of  the  long  period 
motion  observed  even  over  the  shot  point  on  the  high  gain 
refraction  equipment  might  be  considered  to  give  support  to 
a  hypothesis  of  this  nature  for  the  generation  of  Rayleigh 

waves. 

A  deep  hole  drilled  to  38.1  meters  for  velocity  informa¬ 
tion  was  also  used  to  gather  information  on  Rayleigh  wave 
generation.  One  and  one-half  pound  shots  were  fired  at  38.1 
meters  and  at  6.1  meter  intervals  up  the  hole  to  a  depth  of 
7.6  meters.  These  shots  were  recorded  with  the  three  com¬ 
ponent  Sprengnethers  at  2.5#  7.5#  12.5#  25#  37.5#  50  meters 
from  the  hole.  The  P-wave  from  the  shots  has  been  discussed 

previously. 

The  radial  component,  even  at  2.5  meters,  shows  a  pro¬ 
nounced  arrival  shortly  after  the  P  arrival  on  the  vertical, 
too  early  to  be  an  S-wave.  This  event  is  seen  on  Figure  59. 
Its  frequency  (30  cps)  is  lower  than  that  of  the  P-wave 
(59  ops).  At  longer  ranges,  at  which  P  has  a  radial  com¬ 
ponent,  this  event  is  combined  with  the  P-wave.  There  is 
also  motion  on  the  transverse  traces  at  this  time.  The 
nature  of  this  motion  is  not  known  at  this  time . 

A  second  event,  present  on  both  the  vertical  and  radial 
components  but  dominantly  radial  is  also  observed.  Cor¬ 
relation  of  phases  indicates  that  this  is  the  motion  which 
corresponds  to  the  Rayleigh  wave  at  larger  distances.  A 
computer  program  to  calculate  the  travel  tlmeB  of  all 
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direct,  reflected,  and  refracted  events  for  two  layers  over 
a  half  space  was  developed  to  get  accurate  arrival  times. 

These  calculations  indicate  that  this  event  has  the  arrival 
time  of  a  direct  S-wave.  The  identification  of  this  as  a 
direct  S-wave  is  quite  important  because  it  indicates  that 
these  contained  explosions  do  generate  sizeable  SV  motion, 
a  wave  type  that  had  not  been  identified  previously  at  this 
site,  although  they  have  been  observed  at  Alpha  and  Augusta. 
The  maximum  amplitude  of  this  event  at  12.5  meters  shown  in 
the  figure,  is  almost  as  large  as  the  maximum  P  amplitude. 

Figure  60  shows  the  radial  and  transverse  components  of 
ground  motion  over  the  hole  from  0.5  pounds  at  3  meters.  The 
arrival  time  on  both  components  is  that  for  an  S-wave.  The 
motion  is  oscillatory  in  nature  as  it  was  on  the  deeper  shots, 
and  the  rectilinear  character  of  the  particle  motion  may  be 
seen. 

Even  though  this  motion  converges  to  the  Rayleigh  mo¬ 
tion  it  cannot  be  implied  that  it  gives  rise  to  the  Rayleigh 
wave.  The  fact  that  the  Rayleigh  velocity  at  the  frequen¬ 
cies  generated  is  approximately  that  of  the  8-wave  in  the 
upper  layers  at  this  site  would  explain  this  correlation 
without  making  the  assumption  that  the  S-wave  actually  gives 
rise  to  the  Rayleigh  wave. 

At  the  Suffield  test  site,  only  two  shots  of  varying 
depth  were  recorded  on  a  profile.  The  Rayleigh  wave  from  the 
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shallower  shot  had  leas  amplitude  than  from  the  deeper  shot. 
The  P-wave  showed  a  similar  effeot.  This  Is  probably  the 
result  of  a  more  efficient  transfer  of  energy  from  the 
deeper  shot  even  though  the  two  shots  were  of  equal  size. 

It  Is  known  that  the  Rayleigh  wave  from  a  surface  shot 
should  exhibit  a  reversal  of  polarity  with  respect  to  that 
from  a  deeply  burled  explosion.  This  was  found  to  be  the 
case  at  Suf field.  However,  at  Florissant  the  Rayleigh  mo¬ 
tion  from  near-surface  cratering  shots  and  shots  In  water- 
filled  pits  had  the  same  polarity  as  that  from  the  buried 
contained  explosions.  The  reason  for  the  persistence  of 
the  Rayleigh  wave  polarity  with  depth  of  the  source  is  not 
known  at  this  time. 

The  reversal  of  polarity  has  been  clearly  observed  In 
two-dimensional  model  studies  of  exploBlon-generated  waves 
under  Contract  AF  19(628)-l689.  In  those  experiments,  the 
reversal  seems  to  occur  at  a  depth  nearly  equal  to  the 
radius  of  the  zone  of  rupture  produced  by  the  explosion. 
Effect  of  Shot  Point  Medium  on  the  Seismic  Signal. 

The  medium  surrounding  an  explosion  plays  a  dominant 
role  In  determining  the  fraction  of  the  energy  released  by 
the  explosion  that  appears  In  the  form  of  elastic  waves, 
and  the  characteristics  of  these  waves.  The  manner  in 
which  the  medium  controls  the  form  of  the  congressional 
waves  has  been  discussed  In  the  previous  section.  In  that 
discussion,  changes  produced  by  relatively  small  variations 
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in  material  properties  associated  with  different  depths  at 
a  single  test  site  were  considered.  In  this  section  the 
properties  of  the  signal  produced  by  shallow  shots  in  dis¬ 
tinctly  different  media  will  be  examined. 

In  view  of  the  fact  that  the  exploration  seismologist 
has  long  known  that  the  shot  point  material  greatly  affects 
the  seismic  signal,  it  is  surprising  that  so  little  re¬ 
search  on  this  subject  has  been  published.  Adams  and  Swift 
(1961)  have  examined  the  relative  output  of  TNT  charges 
tamped  in  tuff  and  salt,  but  their  observation  that  the 
signal  in  tuff  is  greater  than  that  in  salt  must  be  inter¬ 
preted  in  the  light  of  the  work  of  Nicholls  (1962)  (KIbs- 
linger,  1963c).  Since  the  transmission  of  energy  from  a 
chemical  explosion  is  by  means  of  the  interaction  of  a 
shock  wave  through  the  explosive  material  with  the  bounaary 
of  the  hole  in  which  the  charge  is  placed,  the  ratio  of  im¬ 
pedances  in  the  two  substances  is  very  Important  in  deter¬ 
mining  how  much  energy  is  coupled  into  the  earth.  Results 
obtained  for  one  chemical  explosive  in  several  media  can¬ 
not  be  extended  to  other  explosives  of  greatly  different 
properties  in  the  same  media.  This  point  is  well  covered 
in  a  recent  work  by  Anderson,  et  al  (1963)*  Work  with 
nuclear  explosions,  which  act  as  true  point  sources  in 
that  there  is  no  detonation  wave  and  the  initial  energy 
density  is  very  great,  indicates  that  the  compressional 
wave  amplitudes  are  in  the  ratio  of  Is  0.25s  1.11s  1.6l 
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for  tuffj  alluvlumi  granitei  salt  (Werth  and  Herbet,  1962). 

Assuming  that  the  theory  developed  by  Sharpe  (I942)and 
used  in  the  previous  section  adequately  describes  the  prop¬ 
erties  of  the  P-wave  In  terms  of  the  Inner  radius  of  elas¬ 
tic  response  and  the  elastic  properties  of  the  medium,  one 
must  consider  the  prior  question  of  what  properties  of  the 
medium  determine  the  size  of  this  equivalent  cavity  and  the 
amplitude  of  the  stress  at  Its  boundary  for  a  given  yield. 
This  leads  to  consideration  of  the  very  complex  processes 
occurring  In  the  highly  stressed  region  around  the  Bhot. 

In  view  of  the  lack  of  very  much  published  data  per¬ 
tinent  to  this  question.  It  waB  Judged  worthwhile  to  organize 
the  observations  that  were  collected  during  this  research, 
and  to  present  these  as  observations,  with  no  attempt  at  a 
quantitative  Interpretation. 

The  best  starting  point  is  a  comparison  of  records 
under  similar  circumstances  In  the  four  media  in  which  exper¬ 
iments  were  carried  out.  Tracings  of  the  field  records, not 
reduced  to  true  ground  motion,  are  shown  In  Figure  6l.  The 
static  magnification  for  each  trace  is  given  in  parentheses 
following  the  designation  of  the  component.  Except  for  the 
sandstone  record,  the  distances  are  the  same.  The  depths 
are  very  similar,  and  the  charge  weights  range  from  0.5  to 
2  pounds.  In  general,  the  two  shots  In  unconsolidated 
material  are  similar  to  each  other  in  frequency  content  and 
total  duration  of  motion.  The  records  In  rock  are  charac- 
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F LORIS SANT  -  CLAY  ALPHA  —  LIMESTONE 

.5lb  at  1.3m,  x«30m  1.5 lbs  at  I.Om,  x«30.5m 


SUFFIELD  ALLUVIUM  AUGUSTA  “SANDSTONE 

2  lb*  at  '-8  m,  x  =  30.5  m  I  lb  at  1.2  m,  x=  21.4  m 


FIGURE  61 

SEISMOGRAMS  AT  THE  FOUR  TEST  SITES 
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terized  by  higher  frequencies  and  shorter  duration.  Even 
at  these  short  ranges,  the  records  In  soli  are  dominated  by 
the  surface  wave  motion,  which  is  less  well-developed  for 
the  sandstone,  and  smaller  yet  In  limestone.  At  distances 
of  the  order  of  100  meters,  the  body  waves  from  these  small 
explosions  are  too  small  to  be  detected  with  this  low-gain 
instrumentation.  The  Rayleigh  waves  discussed  in  Chapter 
III  are  recorded  well  in  both  soil  media,  with  dispersion  as 
described  previously.  The  Rayleigh  waves  in  rock  materials 
are  barely  detectable .  In  the  limestone  the  Rayleigh  wave 
appears  as  the  typical  non-dispersed  pulses  expected  for  a 
point  source  near  the  surface  of  a  half-Bpace. 

The  properties  of  the  signals  in  clay,  sandstone,  and 
limestone  have  been  compared  in  some  detail.  Prom  the  ob¬ 
served  velocities  and  densities,  the  moduli  in  these  media 
are  in  the  ratio  of  It  20: 350.  In  the  following  comparisons 
of  response,  rigidity  is  called  upon  as  the  significant 
parameter,  but  other  properties  such  as  yield  stress  or  vis¬ 
coelasticity  may  actually  be  of  equal  importance. 

The  characteristics  of  body  and  surface  waves  observed 
at  Florissant,  Augusta,  and  Alpha  are  tabulated  in  Tables 
5,  6  and  7.  In  this  simplified  analysis  of  the  data,  no 
attempt  was  made  to  separate  geometric  attenuation  from 
absorption  losses  (this  was  done  for  the  Rayleigh  waves  at 
Florissant  and  Suffield,  Chapter  III).  Amplitudes  were 
assumed  to  follow  a  law 
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u  =  U0  x"n 

and  the  approximate  value  of  n  for  each  wave  type  determined. 
This  value  is  listed  under  "Attenuation  Exponent"  in  the 
Tables.  In  general,  the  body  waves  exhibit  an  exponent  be¬ 
tween  2  and  3,  while  the  value  for  surface  waves  is  1.3. 

The  nature  of  the  waves  at  Florissant  and  Suffield 
have  been  discussed  in  the  preceding  chapters,  as  has  the 
SH  motion  at  all  of  the  sites.  Only  a  few  comments  need 
to  be  added  to  the  information  given  in  the  Tables.  At 
Alpha  there  is  no  evidence  of  layering  at  depths  shallow 
enough  to  affect  the  recorded  signals.  The  direct  P,  with 
primarily  radial  motion,  and  direct  S  wave  are  the  only  body 
waves  identified.  The  S-wave  exhibits  both  SV  and  SH  com¬ 
ponents,  in  phase.  The  values  of  corapressional  and  shear 
velocity  give  a  value  of  0.27  for  Poisson's  ratio.  The 
velocity  of  the  non-dispersed  Rayleigh  pulse  is  concordant 

with  the  body  wave  velocities. 

Although  the  sandstone  section  was  thought  to  be  uni¬ 
form  in  composition  to  a  considerable  depth,  a  refracted 
p-wave  was  recorded  from  a  boundary  at  a  depth  of  15  meters 
The  dispersion  of  the  indistinct  Rayleigh  pulse  confirms 
the  presence  of  some  layering.  A  direct  S-wave  was  iden¬ 
tified  at  this  site,  and  a  refracted  S-wave  arrival  was 
also  picked.  The  SH  motion  on  the  record  selected  for  re¬ 
production  in  Figure  61  is  quite  small,  but  at  other  azi- 
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muths,  motion  on  this  component  with  amplitudes  as  large  as 
the  vertical  or  radial  motion  was  recorded. 

The  following  conclusions  are  drawn  from  all  of  the  data 
examined.  They  pertain  to  small  charges  of  dynamite  In  the 
particular  media. 

1)  In  a  dense,  brittle,  rigid  medium,  the  maximum 
ground  motion  Is  smaller  than  In  an  unconsol¬ 
idated  material! 

2)  The  frequency  of  all  events  is  higher  In  the 
rigid  material; 

3)  No  difference  In  the  attenuation  of  surface 
waves  was  found,  while  the  direct  P  wave  is 
attenuated  less  In  the  rigid  material; 

4)  The  ratio  of  P  to  Rayleigh  wave  amplitudes  Is 
greater  In  the  more  rigid  material; 

5)  Shear  waves  (body  waves  and  Love  waves)  are 
bigger  relative  to  P  waves  and  Rayleigh  waves 
in  the  rigid  material.  However,  there  are 
nodal  lines  along  which  the  shear  motion  Is 
very  small; 

6)  The  radiation  of  P  and  Rayleigh  waves  departs 
much  farther  from  symmetry  in  brittle  materials 
than  It  does  in  soft,  compressible  soil. 


OBSERVED  EVENTS  AND  MEASURABLE  CHARACTERISTICS  IN  CLAY 
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APPENDIX 


LARGE  YIELD  EXPLOSIONS  RECORDED  DURING  THE  PROJECT 

Two  large  yield  explosions  were  recorded  during  the 
program,  neither  of  which  met  completely  the  requirements 
to  provide  data  compatible  with  the  seismograms  from  the 
numerous  small  shots.  A  shot  of  100  tons  of  TNT  in  the 
form  of  a  hemisphere  was  detonated  by  the  Canadian  Defense 
Research  Board  on  the  Watching  Hill  Test  Site  at  the  Suf- 
field  Experimental  Station  on  August  3>  1961.  The  experi¬ 
ment  was  conducted  in  an  expert  manner  by  the  Suffleld  staff 
and  excellent  seismograms  were  obtained,  but  because  the 
charge  was  detonated  on  the  free  surface  it  proved  to  be  a 
poor  source  for  the  type  of  studies  that  were  included  in 
this  research. 

The  results  of  this  work  were  reported  to  the  Defense 
Atomic  Support  Agency  and  the  principal  conclusions  were 
given  in  the  Semi-Annual  Technical  Report  No.  2  (Kissllnger, 
1962).  The  entire  set  of  seismograms  has  been  published  as 
Trial  Record  No.  435  by  the  Suffleld  Experimental  Station 
(Jones,  et  al,  196l).  This  shot  did  not  produce  any  unex¬ 
pected  seismological  results.  Several  refracted  P  waveb 
were  identified  that  agreed  in  general  with  the  results  of 
an  independent  conventional  refraction  survey.  Refracted 
SH  waves  from  the  same  horizons  were  recorded,  adding  to 
the  evidence  that  such  events  are  a  normal  effect  of  an  ex- 
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plosion  that  strongly  deforms  the  earth  materials  around 
it.  Analysis  of  surface  waves  was  difficult  because  of  the 
overwhelming  effect  of  the  air  blast  arrival  and  the  ground 
motion  coupled  to  it. 

The  second  large  shot  recorded  was  made  up  of  19,000 
pounds  of  Tovex,  buried  in  a  pattern  of  nine  holes,  30 
meters  apart.  The  shot  point  was  about  26  km  east  of  Pop¬ 
lar  Bluff,  Mo.,  near  the  epicenter  of  an  earthquake  of 
March  3 »  1963*  The  motion  was  measured  on  the  three-compo¬ 
nent  Instruments  at  0.9  km,  and  at  Blx  positions  about  3  km 
from  the  source.  This  shot  was  set  up  and  fired  by  the 
U.  S.  Geological  Survey  on  June  28,  1963  in  order  to  provide 
travel-time  data  for  a  portion  of  the  mid-continent.  There 
was  no  time  before  the  expiration  of  the  contract  to  analyze 
these  seismograms,  but  they  are  of  excellent  quality  and 
will  be  worked  on  In  the  future.  The  pattern  of  charges  was 
set  up  to  diminish  the  close-in  surface  waves,  so  that  it 
was  obviously  not  an  ideal  experiment  from  the  viewpoint  of 
this  study. 

A  well-developed  Rayleigh  wave  was  recorded  at  all 
stations,  as  well  as  a  very  late  prominent  Love  wave.  The 
seismogram  at  the  station  3.2  km  to  the  north  is  of  a 
different  chfiracter  than  all  the  others,  even  the  one  at 
0.9  km  to  the  north,  in  that  it  contains  short  period  Ray¬ 
leigh-type  waves  (probably  higher  mode  events)  that  are 
absent  on  the  other  instruments.  A  conventional  refraction 
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record  along  a  north-south  profile  that  Included  the  sta¬ 
tion  at  0.9  km,  and  extended  almost  to  the  3.2  km  station, 
shows  strong  evidence  of  a  shallow,  burled  fault,  upthrown 
to  the  north,  with  a  throw  of  about  100  meters,  between 
the  source  and  this  north  station  (but  north  of  the  0.9 
km  station).  If  this  fault  really  exists.  It  would  un¬ 
doubtedly  have  a  strong  effect  on  the  short  period  surface 
waves. 

Practically  nothing  Is  known  about  the  tectonics  In 
this  entire  region  of  southeast  Missouri,  a  region  which  Is 
the  source  of  many  earthquakes.  It  Is  Intended  to  follow 
up  this  lead  by  further  geophysical  Investigation,  employ¬ 
ing  the  magnetometer  and  gravity  meter,  as  well  as  further 
seismic  work. 
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